Disordered d-wave superconductors with interactions 



Luca Dell'Anna 1 

1 International School for Advanced Studies (SISSA), 1-34014 Trieste, Italy, 
Dipartimento di Fisica, Universita di Roma "La Sapienza" , 1-00185 Roma, Italy, 
Max-Planck-Institut fur Festkorperforschung, D-70569 Stuttgart, Germany, 
Institut fur Theoretische Physik, Heinrich-Heine- 
Universitat, D-40225, Diisseldorf, Germany 

Abstract 

We study the localization properties of disordered d-wave superconductors by means of the 
fermionic replica trick method. We derive the effective non-linear cr-model describing the diffusive 
modes related to spin transport which we analyze by the Wilson-Polyakov renormalization group. 
A lot of different symmetry classes are considered within the same framework. According to 
the presence or the absence of certain symmetries, we provide a detailed classification for the 
behavior of some physical quantities, like the density of states, the spin and the quasiparticle 
charge conductivities. Following the original Finkel'stein approach, we finally extend the effective 
functional method to include residual quasiparticle interactions, at all orders in the scattering 
amplitudes. We consider both the superconducting and the normal phase, with and without chiral 
symmetry, which occurs in the so called two-sublattice models. 

PACS numbers: 74.20.-z, 74.25.Fy, 71. 23. An, 72.15.Rn 
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I. INTRODUCTION 



After the discovery of superconductivity at high temperatures in doped cuprate materi- 
als Q], a lot of efforts have been devoted to find a theoretical justification of this amazing 
property and to explain the complex phenomenological aspects of these "high T c supercon- 
ductors" . 

One of the novel features of these cuprates is the presence of gapless Landau-Bogoliubov 
quasiparticle excitations in the superconducting phase, due to the d-wave symmetry of the 
order parameter. This peculiarity makes these materials a suitable play ground for studying 
the role of disorder in gapless superconductors. 

Moreover, since cuprates are essentially two-dimensional, quantum interference and local- 



ization should strongly affect the quasipartic 
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The present work, which revisits and extends the analysis in Ref . . would like to 
contribute to the widely studied but somehow still controversial topic of localization effects 
in these unconventional superconductors. There are several theoretical issues which are still 
open and motivate our investigation. 

The first issue concerns the charge transport. Within a BCS description of the super- 
conducting phase, charge is not a conserved quantity due to the gauge symmetry breaking. 
Therefore, although quasiparticles are gapless for a ci-wave order parameter, nevertheless 
long-wavelength quasiparticle charge fluctuations are not diffusive in the presence of disor- 
der, unlike in a normal metal. The diffusive modes just carry spin or energy, which remain 
conserved quantities. It is well known that, as disorder increases, quantum interference may 
lead to the Anderson localization in a normal metal j^, 44, 4^. While it is clear that such 
a phenomenon may suppress spin and thermal conductivities in a <i-wave superconductor, 
the effects on the quasiparticle charge conductivity are not as settled yet. For this reason we 
extend existing quantum field theory approaches built to deal with the truly diffusive modes 
in a ci-wave superconductor to include the charge modes, which acquire a mass term by the 
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onset of superconductivity. In this way we are able to calculate how the charge conductivity 
is modified by disorder in comparison with spin and thermal conductivities. We find the 
same corrections for all of them in all the cases under study but one, when time reversal 
symmetry is broken at half filling, although the conductivities are different at the Born level. 

A second issue has to do with some controversial results about the quasiparticle density 
of states at the chemical potential in the presence of disorder. Within the self consistent 
T- matrix approximation scheme it was found that the density of states vanishes 

linearly as the Fermi energy is approached in the pure system, and acquires instead a finite 
value in the presence of disorder. This was later used as the starting point to build up a stan- 
dard field-theoretical approach based on the non-linear cr-model to cope with the quantum 
interference corrections not included within the self-consistent T- matrix approximation . 
This approach predicts localization and, eventually, a vanishing density of states. 

This standard perturbative technique was nevertheless somewhat unsatisfactory. It was 



argued 2l|, |22j that systems with nodes in the spectrum require a more careful analysis. 



The quasiparticle spectrum in a <i-wave superconductors can be described by 2- 
dimensional (2D) Dirac fermions, with conical spectrum. In the absence of interactions, 
the 2D quantum problem in the presence of disorder becomes effectively a 2D classical prob- 
lem with the retarded-advanced frequency of the single particle Green's function playing 
the role of an external field. This is the non-linear a-model description j^J. On the other 
hand, a classical zero frequency model in 2-dimensions with conical spectrum is analogous 
to a quantum problem of Dirac fermions in 1+1 dimension. Within this scheme the DOS 
is found to vanish (with a different behavior from the approach above) and no localization 
is predicted |22j. In the language of Dirac fermions in 1+1 dimension, the disorder average 
within the replica trick method generates an effective interaction among the one-dimensional 
(ID) fermions, with all the complications that are known to occur. For instance, translated 
in the ID language, the self-consistent T-matrix approach which generates a finite density 
of states at the Fermi energy is analogous to the Hartree-Fock approximation, which al- 
ways leads to density-wave order parameters. However, it is known that Hartree-Fock fails 
completely in ID, which poses serious doubt about the validity of the T-matrix approach 
even as a starting point of a perturbative treatment. From this point of view the contro- 
versy concerns, more deeply, the question of which model correctly describes the quantum 
interference effects. Indeed in the peculiar case in which at most pairs of opposite nodes are 
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coupled by disorder, the perturbation theory around the T-matrix saddle point solution does 
not contain any small parameter, like the inverse conductance in the conventional Anderson 
localization. However, in the most general case of disorder, when all four nodes are coupled, 
we find that a small parameter still exists being related to the anisotropy of the Dirac cones, 
suggesting a conventional field-theory treatment. Looking more carefully at the problem we 
find that the above controversial results are related within the non-linear cr-model to the 
presence (when only opposite nodes are coupled by disorder) or the absence (when all four 
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nodes are coupled) of a Wess-Zumino-Witten term 

A third interesting issue is the role of the nesting property in these kind of systems. Given 
a generic eigenfunction with energy E and amplitude </>(i)e at site i = (n,m), the operator 

0^ii) E = (-l) n+ ™0(i) E , (1) 

which shifts the momentum by (tt, it), generates the eigenfunction with energy —E if nesting 
occurs. This implies an additional symmetry, which has the form of chiral symmetry, at 
E — 0, when the two wavefunctions (1 ± (— l) n+m )(f)E^o+ defined on different sublattices, 
with n + m even or odd, are both eigenvectors. The nesting property occurs when the 
operator O n anticommutes with the Hamiltonian, which is possible in models in which 
the Hamiltonian contains only terms which couple one sublattice with the other, so called 
two-sublattice models. In addition, the chiral symmetry further requires half-filling. Both 
conditions are quite strict and do not represent a common physical situation. Nevertheless 
chiral symmetry leads to quite different and somehow surprising scaling behaviors that are 
worth to be studied. It was seen, for instance, that this symmetry drastically changes the low 
energy density of states. Several models presenting a chiral symmetry were found to have 
isolated delocalized states at the band center at low dimensions. It was argued 0, El that 
these models corresponds to a particular class of non-linear cx-models and it was shown that 
quantum corrections to the (3 function which controls the scaling behavior of conductivity 
vanish at the band center at all orders in the disorder strength, leading to a metallic behavior 
at that value of chemical potential. Moreover, the j3 function of the density of states was 
found to be finite, unlike in the standard Anderson localization. These scaling laws generate 
a divergent behavior at low energy of the density of states. The anomalous terms in the 



action, when chiral symmetry 



rolds, were found to be connected with fluctuations of the 



staggered density of states |52j ]. The modes representing these fluctuations are massive 



in standard non-linear cx-models, while they become diffusive in two-sublattice cases. For 
this reason in the conductance channels with both positive and negative frequencies acquire 
diffusive poles and contribute to quantum interferences corrections. This is what we find also 
in our two sublattices ci-wave superconductive model that presents extended states at the 
band center which are associated with diffusive spin transport. The DOS diverges and shows 
a behavior similar to that of 5^]. Furthermore, we find an unexpected charge conductance 
behavior. As we have said before, although charge modes in ci-wave superconductors are not 
diffusive, nevertheless quantum interference corrections affect the charge conductance. In 
particular, when chiral symmetry holds but time-reversal symmetry is broken, quasiparticle 
charge conductivity is suppressed, but spin and thermal conductivities stay finite, leading 
to a spin-metal but charge-insulator quasiparticle behavior. Moreover we saw that, even 
though magnetic fields or magnetic impurities introduce on-site terms in the Hamiltonian 
that spoil the full sublattice symmetry, staggered fluctuations are not totally suppressed 
introducing other symmetries in the model under study. We saw, for instance, that the 
problem of ci-wave superconductors with chiral symmetry and magnetic impurities can be 
mapped to a U(2n) non-linear cx-model and belongs accidentally to the same universality 
class as ci-wave superconductors far from the nesting point embedded in a constant magnetic 
field. 

From the point of view of the cuprate ci-wave superconductors, it is not unlikely that 
chiral symmetry may play some role, especially in underdoped systems close to the half- 
filled Mott insulator. Indeed it is believed that the impurity potential is close to the unitary 
scattering limit |2fJ which, by taking out one site, reduces to a random nearest-neighbor 
hopping. Furthermore, although the band structure does not have a perfect nesting, the 
superexchange interaction which stabilizes a Neel antiferromagnetic phase at half-filling may 
effectively reduce the energy scale at which deviations from perfect nesting get appreciable. 

The last issue that we will discuss is the role of the residual quasiparticle interaction and 
its effects on the conductivity and on the density of states. Following the original Finkel'stein 
approach j^J, which extended the effective functional method to the disordered electron- 
electron interacting systems, we introduce the effective quasiparticle scattering amplitudes 
in different channels, firstly considering systems without sublattice symmetry. We find 
that, consistently with the charge not being a conserved quantity, the singlet particle-hole 
channel does not contribute while the triplet channel does. On the other hand, the scattering 
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amplitude in the Cooper particle-particle channel acquires a factor 1/2 with respect to the 
normal metal state, which corresponds to the fact that only the phase of the order parameter 
is massless. We see that the effective interaction, which we assumed being repulsive, has 
a delocalizing effect enhancing the density of states. We extend the Finkel'stein model in 
order to include the nesting property, by introducing additional scattering amplitudes with 
(tt, 7r) momentum transferred, and we evaluate the new corrections to the density of states 
and to the conductivity. We also consider the metallic phase which can be of relevance for 
2D metals or semimetals. Moreover, we notice an interesting fact which occurs at half-filling 
for a two-sublattice model in both the superconducting and normal phase. Depending on 
the sign of the interaction the staggered particle-hole fluctuations, being diffusive, can lead 
to a log-divergent staggered susceptibility, and a Stoner instability towards a spin or charge 
density wave. 



II. THE MODEL 



The characteristic feature of a (i-wave superconductor is the existence of four nodal points 
where the order parameter vanishes. To study the low temperature transport properties of 
a (i-wave superconductor, we consider the following model defined on a two-dimensional 
square lattice with lattice constant a: 

n = - E E (^ c *W + H - c ) + E ( c It4 + c ]A) + H - c ] > ( 2 ) 

(ij) <* (ij) 

where (ij) means that the sum is restricted to nearest neighbor sites, c ia creates an electron 
with spin a =|, I at site i, while c ia annihilates it. We take a gap function A;.,- of (i-wave 
symmetry. The hopping matrix elements are independent random gaussian variables with 
average value t and variance ut, and satisfy £^ = tji e R, and (fry = —<pji, with 0^ zero 
or finite depending whether or not time reversal invariance holds. The spectrum of the 
Hamiltonian (J2J) possesses a nesting property. This implies an additional symmetry (chiral 
symmetry) at half filling [^. The localization properties are quite different whether chiral 
symmetry holds, which corresponds to the Fermi energy Ep = (half-filling), or broken, 
Ep 7^ 0. In the latter case, the localization properties of (j2J) are analogous to models in which 
on-site disorder is present or next-nearest neighbor hopping is included, which break chiral 
symmetry everywhere in the spectrum. For this reason while dealing with a bipartite lattice 



that induces an higher degree of symmetry we can reduce the problem to the standard 
ci-wave case only by introducing an on site term in the Hamiltonian, thus spoiling chiral 
symmetry. 

In the absence of randomness the quasiparticle spectrum has four nodes at (±kp, ±kp, )■ 
In the vicinity of each gap node the Fourier transform of —tji, namely e& = — 2t cos(k x a) — 
2tcos(k y a), varies linearly perpendicularly to the Fermi surface while the Fourier transform 
of Ay, that is Afc = 2 A (cos(k x a) — cos(k y a)), varies linearly parallel to the Fermi surface. 
Let us rotate the axes from k x ,k y to ki,k 2 by 7r/4 rotation, and define a Fermi velocity, 
Vi perpendicular to the Fermi surface, and a gap velocity v 2 parallel to the same surface. 
Then, close to the nodes the quasiparticle spectrum is 

E k ~ yjv*k* + v*kl (3) 

for nodes along k\ axis and 

E k ~ sjvlkl + v\kl (4) 

for nodes along k 2 axis. The spectrum, in the vicinity of each gap node takes the form 
of a Dirac cone whose anisotropy is measured by the ratio of the two velocities. As we 
will see afterward a strong anisotropy brings in a weak coupling regime the non-linear a- 
model representative of our disordered system, making the perturbation theory and the RG 
approach suitable tools of investigation. 

We now analyze the disordered Hamiltonian Q by using the replica trick method within 
the path integral formalism 3|. We introduce the vector Grassmann variables q and q 
with components Ci )0VPj0 and c^a,p,ai where i refers to a lattice site, a to the spin, p = ± is 
the index of positive and negative (—u)) frequency components, and a = 1, ... ,n is 

the replica index, as well as the Nambu spinors 

^ = 4= I Ci 1 and #,= [C^]*, (5) 
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with the charge conjugation matrix C = %a y T\. Here and in the following, the Pauli ma- 
trices a;, (h = x,y,z) act on the spin components, s& {b = 1,2,3) on the frequency (re- 
tarded/advanced) components, and r& (b = 1, 2, 3) on the Nambu components c and c. The 
action corresponding to (j2J is 

S = ( '°' i7 i + iA^-raSi - i^wsg) ^ (6) 



where the source term uj^/iS^i is introduced in order to reproduce positive and negative 

frequency propagators. As in the standard Abrikosov-Gorkov-Dzyalozinskii approach 

to superconductivity, the gap function couples fermions with opposite frequency uj. 

If magnetic impurities are present, we must add to Eq. (jUJ) an additional spin-flip scattering 

term 

-^^fiUiTza- S^i, (7) 

i 

being Si the impurity spin and Ui the corresponding random potential. The same term, with 
Si = B and = B, gives the Zeeman splitting in the presence of a constant magnetic field 
B, which also breaks time-reversal invariance. 

If an on site term is present or if we are far from half filling, we must add 

-J^Mi, (8) 

i 

which spoils chiral symmetry, fii being a random variable or Ep respectively. 

III. DISORDER AVERAGE 

We derive the effective quantum field theory for the disordered <i-wave superconduct- 
ing model described in the previous section, following the work by Efetov, Larkin, and 
Khmel'nitsky 5J|. As usually, this field theory is a non-linear cr-model where the broken 
gauge symmetry enters as a reduction of the symmetry of the Q-matrix fields with respect 
to a normal metal. 

The starting generating function is 

Z = I V^VmVtP(t) e~ s , (9) 



where P(t) is the gaussian probability distribution of the random bonds. Within the replica 
trick method, we can average e over disorder, obtaining 

Z = I V^Vme~ s ^i. (10) 



The effective action S e ff — So + Si mp is the sum of a regular part So given by Eq.(jnj) with 
tij = t (and (pij = 0) plus the impurity contribution 

S imp = -2uh 2 (^*i) 2 = -2«V Yl (*i*0 ( n ) 

(y) (ij) 
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since fy^j = ^j^i- By introducing the fields 

Xf = *?*f, (12) 

where a and (3 is a multilabel for Nambu, advanced-retarded and replica components, we 
can write 

S mp = 2u 2 t 2 XfX^ a = 2uH 2 J2 Tr (X^) . (13) 
In Fourier components it becomes 

S imp =^ J2w g Tr(X q X_ q ), (14) 

qeBZ 

where BZ means the Brillouin zone, V the volume and 

w g = 2u 2 t 2 (cos q x a + cos q y a) , (15) 

a being the lattice spacing. We can decouple Eq. (JHJ) by an Hubbard-Stratonovich trans- 
formation, introducing an auxiliary field. However, since w q = — w g+ ( 7rj7r ) and w q > if q is 



restricted to the magnetic Brillouin zone (MBZ), we neec 
defined within the MBZ, Q 0q = Q\_ q and Q 3q = Q\_ q j^l , 



;o introduce two auxiliary fields 



521 ] . through which 



Simp — T7 ~ Tr [QoqQo-q + Q3qQ?>-q] 

V q eMBZ 4w « 

- % y E Tr [Qo^~ q + iQs q X- q -(^)} ■ ( 16 ) 

qGMBZ 

The above expression shows that Qo corresponds to smooth fluctuations of the auxiliary 
field, while Q3 to staggered fluctuations. Namely, in the long- wavelength limit, the auxiliary 
field in real space is 

(), Q 0j ■ ii-iychj- (ir) 

where j is the site. However, to make sublattice symmetry more evident, it is convenient to 
use a unit cell which contains two sites, one for each sublattice. Indicating with R a new 
unit cell vector and with A and B the labels for the two sublattices, we introduce a two 
component field 

*R = I I , (18) 



BR 
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through which we can rewrite Eq. (|17|) in the following way 

Qr = QoRlo + iQ3R73, (19) 

where 75 (b = 1,2,3) are Pauli matrices acting on the vector (|18j) . Notice that, differently 
from the non-chiral case, Q here is not hermitian, in fact 

Qr = QorIo ~ iQmlz = IiQrIi = 72<2i?72, (20) 
since Qo and Q3 are both hermitian. 

IV. SYMMETRIES 

From Eq ()16|) we can deduce that a unitary transformation acting on the spinor — > 
T^ji results into a rotation of the matrix field 

Q -> CT'C'QT. (21) 



Therefore, as originally done by Wegner 46], we analyze the symmetries underling the the- 
ory in order to distinguish the soft or transverse modes from the massive or longitudinal 
modes, once we will define the vacuum states for the Q matrices by the saddle point ap- 
proximation. Since the Hamiltonian parameters couple sites of the two different sublattices, 
we can consider generally two different global unitary transformations, one for sublattice A 
and another for sublattice B 

y A = T A ^ A , m B = T B ^ B . 
For frequency to = 0, the action is invariant if 

CT^^HabTb = H AB , (22) 
being C is the charge conjugation matrix i<y y Ti, that implies 

CTXC'Tb = 1, 

CT B C l T A = 1, (23) 
valid even for non-superconducting states, as well as 

CT^C^Tb = r 2Sl , 

CT t B C t r i 8 1 T A = r 2 si, (24) 
10 



in the presence of a real superconducting order parameter. If time reversal symmetry is 
broken, namely if the hopping parameter acquires a phase e* 9 ^ 7 " 3 , we must further impose 
that 

CT A Crr 3 T B = 7-3, 

CT t B C t r 3 T A = r 3 . (25) 

In the presence of a constant magnetic field B, which introduces a Zeeman term B z r 3 a z in 
the Hamiltonian, we add the following constrains, 

CT 1 a C 1 t 3 T b = r 3 , CT t A C t T 3 a z T A = r 3 a z , 

CT B C t r 3 T A = r 3 , CT t B C t T Z a t T B =T Z <T z , (26) 

and finally 

CTACV3T4 = r 3 , CTACV30T4 = r 3 a, 

CT B C t T 3 T B = T 3 , CT t B C t T Z 3T B = T Z 3, (27) 

— * 

in the presence of magnetic impurities, represented by the term S ■ ar 3 in the Hamiltonian 
with S a random vector variable. 

If we are not at half filling or there are on-site impurities, namely if sublattice symmetry 
does not hold, T A and T B have to satisfy 

CTX&Ta = 1, CT B C l T B = 1, (28) 

in order to be symmetry transformations. Finally, a finite frequency u, which acts as a 
symmetry breaking field, leads to the additional constrains 

CTlC%T A = S3 , CT t B C t s 3 T B = s 3 . (29) 

All the above conditions imply that the unitary matrices T belong to a group G if u = 

which is lowered to a group H at u 7^ 0. 

The unitary transformations, T A and T B , can be written as 

w + w 3 w -w 3 

T A = exp , T B = exp , (30) 



with antihermitian W's. Charge conjugacy invariance, through Eqs. (|23j) . implies 

CWlC 1 = -W , (31) 
CWlC 1 = W 3 , (32) 
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while the presence of superconducting term leads to the following constrains 



[T2Si,W m ] =0. (33) 

If time reversal invariance is broken 

[r 3 ,Wo( 3 )] =0. (34) 
If magnetic impurities are present, then 

[r 3 a,W ] = {r 3 a,W 3 } = 0, (35) 

while 

\ 3 a-B,W 3 }=0, (36) 



r 3 a ■ B, W 

in the presence of a constant magnetic field. Finally, if chiral symmetry is broken by terms 
which couple the same sublattice, we must set W 3 = 0. In the presence of finite frequency 
uj we have further to impose that 

[W ,s 3 } = {W 3 ,s 3 } = 0. (37) 

V. SADDLE POINT 

The full action, after the Hubbard-Stratonovich decoupling, 

S = ~ E f^^o - H° k 5 q0 + 1q_ 9 ) v k+q + i J2 ^ Tr [Q q Q\] » (38) 

k,q ^ ' q q 

by integrating over the Nambu spinors, transforms into 

where if is the regular part of the Hamiltonian. In momentum space it reads 

H k = e k + i A k r 2Sl = E k e 2i6k ™ , (40) 

where 



E k = V e i + A l ( 41 ) 

and 

cos20 fc = ij_, sin2£ fc = ^. (42) 
12 



Let us look for a saddle point Q sp oc o"o which has a tqs^ component £ as well as a T2S1 
component F, both /c independent. Therefore 

G^ 1 = zo;s 3 — e k — zA fc r 2 si + i£s 3 + iFr 2 s l} (43) 

where we introduce explicitly the symmetry breaking term, namely US3. We notice that the 
new pairing order parameter is — F, so that, by defining 



E k = J el + (A k - F) 2 , (44) 



as well as a modified 9 k , we find the self-consistency equations 

u 2 t 2 \ - 
£ = i-g- 2^ Tr (G fc s 3 ) , 

k 

F = z^^Tr(G fc r 2Sl ), 



where 



Therefore, 



Q — Q-iOkViai \ Q -idkT2Si 

-E k + i(u + £)s 3 



F 



(S + ")E ^ 2 , ,1 i ,, 2 > ( 45 ) 
^ V = ^F V , (46) 



2 ^ E 2 + ^ 2 V^ + S2 



where the last identity holds for <i-wave order parameter. Notice that, for s-wave symmetry, 
these equations coincide with those found by Abrikosov, Gorkov and Dzyalozinskii j^]. The 
first equation implies that 

g UH 2 1_ 



which, inserted in the equation for F, leads to 



0. (48) 



S + cu 

Being u; non zero, although infinitesimally small, this equation has the solution F — 0. 
Therefore, only S 7^ and such that 

u 2 t 2 1 
13 



The above self-consistency equation leads to the Born approximation result 

S = tiuH 2 v = jw v oc e — W« 2 * 2 ; (50) 

with v being the density of states at the chemical potential. Finally one gets the saddle 
point value 

Qsp = o"o To s 3 S. (51) 



VI. TRANSVERSE MODES 

Now we will consider the transformations that leave the total Hamiltonian unchanged 
while rotating the saddle point, that is to say the transformations that allows to move from 
a vacuum state to another one. The degrees of freedom of these transformations are the 
Goldstone modes which are massless and whose number is equal to 

dim(G/H) = dim(G) - dim(H), 

where G is the original symmetry group and H is the symmetry group that preserves the 
vacuum. The coset G/H tells us how many generators are broken. Since the saddle point has 
the same algebraic form of the frequency term in the action, the cosets related to Goldstone 
modes are obtained exactly by that transformations which are excluded in the lowering of 
symmetries due to finite frequency. In the following we will denote by T only those kind of 
transformations, represented in sublattice space by 

T = e ^ , (52) 

where the 7's are 2x2 matrices acting on vectors (jlSjl . We adopt this exponential form, 
which is the adjoint representation of the group whose transverse modes belong to, because in 
this way small fluctuations from the saddle point can be easily written by Taylor expansion. 
For each Wq^ we can separate the singlet term from the triplet one, writing 

W = W s + iff ■ W T , (53) 

where the Pauli matrices a a , a = x,y,z, act on spin space. In addition we rewrite W in 
r-components (Nambu space) 

3 

W s = Wsoro + iY^WsjTj, (54) 
3=1 

14 



3 

W T = W T0 T + i ^ W T jTj. (55) 
i=i 

Moreover, for each r-component, we write (a = 0, 1, 2, 3), 

3 

Ws(T)a = W^S(T) (56) 

Each component of in Eq. (JSTIj) is a n x n matrix in replica space. We denote by S 
and A the symmetric and the antisymmetric matrices in replica space and by R and I the 
real and imaginary ones. These matrices are collected in Appendix A, together with the 
corresponding cosets. The transverse components of the IV-fields have to fulfill 

{W ,s 3 }=[W 3 ,s 3 ] = 0. (57) 

VII. NON-LINEAR a-MODEL 

Here we derive the effective field theory describing the long wavelength transverse fluctu- 
ations of Q(R) around the saddle point. In general terms we may parametrize the Q- matrix 
as follows (cfr Eq (J2TJ)) 

Q P {R) = f(Ry [Q sp + P{R)) T(R) = Q(R) + f(RyP(R)T(R), (58) 

where T involves transverse massless fluctuations and P longitudinal massive ones, Q sp = 
S S3 being the saddle point. In Eq (}58|) . we used the short notation 

ft = CT'C 1 = 7l Tt 7l = 72 TV (59) 

Since only the T's are diffusive, at the moment we concentrate just on them, neglecting the 
P's and writing the action in terms of Q(R) = T{R)^Q sp T{R) alone, so that QQ' = Q sp 
is invariant, eventhough a term involving massless modes from integration over massive 
ones might appear. Afterwards we will reconsider this point. By integrating (jHSj) over the 
Grassmann variables, we obtain the following action of Q: 

- S[Q] = ~ J2 ^r Tr + \ Tr ^ - # (0) + iQ] ■ (60) 

<? q 

Since global transformations leave invariant the action with u = 0, the effective theory in 
terms of Q-fie\ds can be found by expanding to the second order in the gradients and to the 
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first order in u. To this end we rewrite the second term of S[Q] as 

^Tr In (iuofs z T ] - TH {Q) T ] + iQ sp ^ . (61) 

H RR , involves either 71 and 72, while T involves 70 and 73. Then 

f(R)H° RRI T(R>y = H RR , + (r(#) f - T(J2)t) H° RR , (62) 
~ H° RRI - f(R)ff(Ry ■ (fi - #) H° RRI 

+ ij^f (R)d ij f(R.y(R i - ^)(^ - i?;.)^, = fl^w + 

with dij = gr^aR. ; -Ri being the components of i?. Differently from the tight-binding Hamil- 
tonian case [52], in which the term (R — R')H RRI is related to the charge current vertex, in 
BCS Hamiltonian, since the charge is not a conserved quantity, the spin current vertex is 
involved, since the spin is conserved. This can be seen writing the continuity equation 

iVJ(R) = [H°, Pspin (R)}, (63) 

with p sp i n (R) = cL-Crj — cLcrj, from which we obtain the following expression for spin 
current on the basis (jHJ) 

J(R) = -1 J2(R ~ Ri)*rH RRi <t z * Ri . (64) 
Ri 

We have chosen z as the spin quantization direction. Using Eq. (|62|) the term ()61|) in the 
action can be written as 



^Tr In (iufs-zT ] - U - H {0) + iQ s 
— Tr In G + ^Tr hi (l + Giufs 3 T^ -GU^j , (65) 



where G = (—H^ + iQ sp ) 1 is the saddle point Green's function. By expanding in lo the 
following term is found 



LO 

i— 
4 



Tr (G TstA = ^-Tr (s 3 Q) . (66) 

V / ZWn 



■'0 

The second order expansion in U contains the terms: 

- l -Tr{GU) (67) 

and 

-^Tr(GUGU). (68) 
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Neglecting boundary terms coming from first derivatives and keeping in ()67|) the component 
of U containing second derivatives, we get 

- \tt [f{R)d l3 f{R)~ l (R t - JO (Rj - R'j) H^G(R', R)} . (69) 

Now let us consider the correlation function 

x,,i(R,R';t,ti,t 2 ) = (t [iW^^^c^Wc^tOJ^^K)^^)]), (70) 

where \i = 0, 1, 2, and 

the spin density matrix elements, while J 1 , for i = 1,2, are the matrix element components 
of the spin current and read 

J Rl R 2 (R) = ~<Ri ~ R2)H° RlR2 a z 5 RR2 . (71) 
By means of the continuity equation in the hydrodynamic limit we obtain the Ward identity 

xoM R 'i E ) = E - R ^ - R 'i) Tr ^ e ) r I',r) ■ ( 72 ) 

RR' RR> 

Through Eq ([72]). Eq. (|p^j) turns out to be 

-I — u 

X 



Trlf^^J?)- 1 }, (73) 



16 

which, integrating by part, is also equal to 

y++ r~ 
- -f^-Tr ^T^diTiR^T^djTiR)- 1 

= —X^TriDtDj). (74) 
Here we have introduced the matrix D(R) with the i-th component 

Di(R) = D ,(R) l0 + D 3>i (R) l3 = f{R)d i f{R)- 1 . (75) 
The second term of the expansion in U, given in Eq. ([68)1. is 

-\tt{gugu) = 1 -Y j Y. Tr {v( R )- a * G ( k ) 3 ( R ) ■ a * G ( k ) } ■ ( 76 ) 

k R 
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Since the Fourier components of the Green's function and of the spin current operator in 
the long- wavelength limit, supposing Uj = tji G R and Ay = A^ G R, can be written in the 
following way 

= . (77) 

Jk = V fe e fc 7 1 cr z + iV fc A fc r 2 Si7i(j 2; , (78) 

and the positive and negative frequency Green's functions are related by G + = siG~S\, we 
have 

J % k o z G + Dj = ^{D j + llS3 D 3 s 3ll )J l k a z G + + ^(Dj - 71*3^371)4^" (79) 
Jl being a component of Jk- Taking advantage of this result, Eq. (fTBj) becomes 

^4 + Tr [A-Dj + DiS^DjS^} + ^xtfTr [B i B i - D lS3ll D r s 3ll ] , (80) 
which, summed to Eq. (fTIj). gives 

^ (XS _ - 4- + )Tr[A^ - A S3 7i^«37i] = -^rrC^gt), ( 8 1) 

where 

^ = -^E Tr [4 (^W - J i ( G+ W ~ G ~W)} ( 82 ) 
k 

is the spin conductivity since Jk is the spin current vertex. For i = j, using Eqs. (|77I78|) . 
we obtain the following quantity 



" TlV 4^ 



Ve fe • Ve fc + VA fc • VA A 



,2 



7T fi^2 



where, as previously defined, t>i and t>2 are the velocities perpendicular and parallel to the 
Fermi surface. This is the quasiparticle conductivity in the Drude approximation 3, [3 
which corresponds to the spin and the thermal conductivities since both energy and spin are 
conserved quantities, therefore fluctuations of energy and spin densities are diffusive. One 
can expect that in the limit in which pair order parameter goes to zero a corresponds also 
to the charge quasiparticle conductivity in Born approximation since in that case the charge 
turns to be a conserved quantity as well. This is what we will find rigorously in Section IX. 
Inserting Eq. (J53)l in Eq. (JHTj) and adding Eq. (JMj) . we arrive finally at the following 
non- linear cr-model action 

S \Q} = ^tW 7 / dRTr(d fM Q(R)a^d u Q^(R))-—Tr(s 3 Q(R)). (84) 

6ZL, J Wq 
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A particular metric a^ v appears, where fi, v = 1,2 denotes the directions k\ and k 2 , de- 
pending on which nodal contribution we include in Eq. (j84j) through Eq. (}8*2*j) . Specifically 



a: 



5^ v for 4 nodes which is the model we have considered starting from Eq. (j2J). If 



we assume to suitably modify the disorder such that it couples quasiparticles belonging to 
only one node or two opposite nodes, namely if we consider only forward scattering in the 
presence of a random potential or including also backscattering processes, the metric to be 
considered is a ul/ = S uu h -iY •> for one node or a uu = S uu > 2 for two opposite nodes. 

At this point it is important to discuss the differences which occur when disorder couples 
at most two opposite nodes, or the most generic case where all nodes are coupled together. 
We anticipate that the logarithmic terms, which appear upon integrating the gaussian prop- 
agator, derive from the expression 

1 f d 2 k 1 -i / \ 

= <7ln(...), (85) 



7TO- J An 2 k^a^ky 

where the effective coupling constant controlling the perturbative expansion is gives by 

g = — — , for 4 nodes, (86) 

2tt 2 (x 

g = I v l±A i f or i no de. (87) 

We readily see that, up to terms of order u A , the disorder strength, for one node g = 1 or 
two opposite nodes g = 0.5, so that the system is never in a weak coupling regime and the 
saddle point physics with weak perturbative quantum interference effects is never realized 



22 1 in which 



not even as a crossover transient. This is the situation analyzed in Ref. 
the authors found a vanishing density of states approaching the chemical potential with a 
disorder strength dependent exponent in the one node case while a universal exponent is 
found when two opposite nodes are coupled by the disorder. In Ref. j^, 57, 58, 59, 0] 
it is pointed out that this results depend on the presence of a Wess-Zumino-Witten term. 
On the other hand, assuming for cuprates the generic 4-nodes case, since experimentally 
t> 2 ~ fi/15, g << 1 and a perturbative expansion in g is still meaningful. In the following 
we will consider the latter situation. 

The full expression of the Q-matrix is expressed by Eq. (|58J) where the massive modes 



are 



P = (Poo-so + ito) 7o + i (P31S1 + P32S2) 7s, (88) 
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being all P's a hermitian. Charge conjugation implies that CP t C t = P. Writing the free 
action of Qp(R) and expanding w q we obtain two contributions, a pure massive term and 
a term where massive and massless modes are mixed. Integrating over massive modes as 
in Ref.0 we find another term which we have to add to the action for the transverse 
fluctuations 

7T 



n / dRTr 



s . 3,^ -- ; ■ • L g T ( J R)Vg( J R) 73 j Tr \Ql(R)VQ{Rft 3 

II is a parameter related to staggered density of states fluctuations 



(89) 



This term exists 



only if chiral symmetry is present. Indeed a further term could appear, namely the Wess- 
Zumino-Witten term which is calculated in detail in Ref. 0| within 

the same approach. It is found that this term accidentally cancels out thanks to the four- 
fold symmetry of the Dirac nodes. Notice that if such a symmetry of the spectrum is 
broken the WZW would appear again and we would take it into account obtaining very 
different scalin g be haviors. However, normally, the system flows to the strong coupling 



regime 
canceled. 



56 



57 



58 



59 



601 ]. where all the nodes are locked and again the WZW term is 



VIII. RENORMALIZATION GROUP 



We now study the scaling behavior of our action by means of the Wilson-Polyakov renor- 
malization group 0, 0]. We will also show how to evaluate the one loop correction to 
the conductivity, namely to the stiffness parameter of modes which acquire a mass, like the 
charge fluctuation inside the superconducting broken symmetry phase or the spin modes 
when spin isotropy is broken. According to the previous section the final expression of the 



action describing the transverse massless modes in the long-wavelength limit is 
S[Q] = ^J dRTr (VQ(R) ■ VQ(i2) f ) - ^- J dRTr (s 3 Q(R)) 



41] 



7T 



8-32 



II / dRTr 



Q\R)VQ{R) l3 Tr Q\R)VQ{R) l3 



(90) 



where we have rescaled the Q matrices by a factor E so that the new Q is normalized to 
unity, i.e. Q^Q = I and Q sp = s 3 . Since Q = T^ s 3 T = s 3 T 2 = s 3 e w , at the gaussian level 
the first term in the action is simply 



7TO" 



32 



dR Tr ( VQ f VQ 



7T<7 

32 



dR Tr {VWVW 



(91) 



20 



while the last term is 



ttII 



32 ■ 8 



dRTr 



Tr 



Q f VQ73 



7TII , , , 

~ / dRTr 

64 



VW, 



Tr 



(92) 



Depending on whether the W components defined by Eqs. (|53H56|) are real or imaginary 
matrices in replica space, which may be either symmetric or antisymmetric (see Appendix 
A), we find the following gaussian propagators for the diffusive modes 

n 



(WSsia(k)Wf Sia (-k)) = ±D(k) (6 ac 6 bd ±6 ad 6 bc ) + D(k)- 



■ 8q3 S ss 5 i0 S a0 S ab 8 cd} (93) 



a + Un 

where q = 0, 3 refers to the homogeneous or the staggered component, S = S,T for singlet 
or triplet term, (J53J) . % — 0, 1, 2, 3 refers to the r-components, (|51 jl .(|55 |l . a = 0, 1, 2, 3 refers 
to the s-components, ()56|) . a, b, c, d are replica indexes, the ± sign in front refers to real (R) 
and imaginary (I) matrices, while the ± sign inside the brackets refers to symmetric (S) or 
antisymmetric (A) matrices, n is the number of replicas and finally 

1 1 



D(k) 



2na k 2 ' 



(94) 



This propagator in two- dimensions will induce logarithmic singularities within the pertur- 
bative expansion. A standard way to handle those divergences is provided by the Renor- 
malization Group (RG). In particular we here apply the Wilson-Polyakov RG procedure 



El, |62j,|63j, which is particularly suitable to handle with the non-linear constraint QQ^ = I. 
By this approach one assumes 

T(R) = Tf(R)T s (R), 

where Tj involves fast modes with momentum q £ [A/s, A], while T s involves slow modes 
with q £ [0, A/s], being A the high momentum cut-off, and s > 1 the rescaling factor (in 
terms of the elastic scattering time, r = ^g, A is given by §;A 2 ~ 1, ^- being the diffusion 
constant at the Born level). The following equalities hold 



Tr 
+2Tr 
+4Tr 



VQtVQ 



Tr 



D sll Q s I) s Q) lx 



D s Q}VQ f 



- 2Tr 



D S D S 



(95) 
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where Qf — TtQ sp Tf and D s = T s VTj, as well as 



Tr 



Tr 



Tr 



VW. + VWA^ -Tr (VWs + VW/Wg 



(96) 



Since the fast and slow modes live in disconnected regions of momentum space, only the 
stiffness Tr VQ^VQ generates corrections. By expanding in Eq. (|95|) the terms which 
couple slow and fast modes up to second order in Wf, after averaging over the fast modes 
within a one loop expansion, the stiffness generates the following contribution to the action 
for the slow modes 



™ / dRTr 
32 



If 



/• 



(97) 



where 



Si 



32 



dRITr 



D-fiQ sp DQ sp Wfri 



-2Tr 



spw f 

D^Q^WfDWfQs^i 



and 



S-2 



32 



dRTr 



DWfVW; 



(98) 



(99) 



After the effective action for the slow modes have been obtained, the slow modes momenta 
are then rescaled back according to : 



q e 



0. 



A 



s 



where q' G [0, A] runs over the original momentum space. In this way the model is mapped 
onto another model defined onto the same range of momenta with renormalized parameters 
a(s) and n(s). The logarithmic singularities are controlled by the following dimensionless 
coupling constants 

9 = 7^, C =A^- (100) 
When chiral symmetry holds, namely when W 3 is massless, the new coupling c has to be 



included. However, the combination a + nil can be shown [50|, |52j to represent the stiffness 
parameter of an abelian degree of freedom connected to the Tr(W 3 ), which is finite and 
commutes with all other degrees of freedom. This implies that o + nil is a constant of the 
RG flow, namely that 

& = — -■ ( 101 ) 

g 2 n 
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Since the theory is well behaved in the n — > zero replica limit, this indirectly proves that 



lim f3„ = 0, 

n— >0 

namely that when chiral symmetry holds and when Tr(W 3 ) is massless the model stays 
metallic with a finite conductance. 

The final results of the RG are collected in the Table in which the (3 functions for g 
(f3 g = dg/dlns) and for the density of states (DOS) (f3 p = din p/ dins) are listed for the 
different symmetry classes. The density of states scaling behavior is obtained through its 
expression in the Q matrix language 

p = ^Tr(s 3 Q), (102) 

which allows a very simple loop expansion. 

In Table |U we also list the coset spaces G/H for the different classes (i) time reversal 

n n n 

time 

reversal symmetry is broken by introducing random phase with ehiraJ symmetry }5l j^lj ] 
or without |32l; (Hi) a magnetic field is applied in the presence of chiral symmetry [4l| or 



without it k3 




651 ]: and fina 



symmetry 4l| or in its absence 



ly (iv) in the presence of magnetic impurities with chiral 



TABLE I: Coset spaces and (3 functions for the coupling g and for the DOS p in the different 



universality classes. V 



c+n g 



and T is the time reversal invariance. 





Coset space 


Pa 


(3 P 


Yes chiral, Yes T 


U(4n)xU(4n)/U(4n) 


8ng 2 


(r/4 - 8n)g 


Yes chiral, No T 


U(4n)/0(4n) 


Ang 2 


(-l + T/4-4n)g 


Yes chiral, magnetic field 


0(4n)/0(2n)xO(2n) 


(2n - l)g 2 


-2ng 


Yes chiral, spin flip 


U(2n)/U(n)xU(n) 


2 

rig 


-ng 


No chiral, Yes f 


Sp(2n)xSp(2n)/Sp(2n) 


2(2n + l)g 2 


(-l-4n)g 


No chiral, No T 


Sp(2n)/U(2n) 


(2n + l)g 2 


(-l-2n)g 


No chiral, magnetic field 


U(2n)/U(n)xU(n) 


2 

rig 


-ng 


No chiral, spin flip 


0(2n)/U(n) 


(n ~ l)g 2 /2 


(1 - n)g/2 



According to Table HJ in the zero replica limit we obtain that, if chiral symmetry is absent 
and for non magnetic impurities, the conductance vanishes, and the DOS, which is finite 
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within the simplest Born approximation, is suppressed. As shown by Ref. [32$, in the 
localized phase the DOS vanishes as \E\ or E 2 depending whether time reversal symmetry 
holds or not. Quite surprisingly, magnetic impurities give a derealization correction to the 
conductance, as well as a DOS enhancement. On the contrary, if chiral symmetry is present, 
the conductance stays finite, or even increases in the presence of a magnetic field. Without 
magnetic field and in the absence of spin flip scattering, the DOS, according to the above 
/3-function, diverges approximately like, p(E) ~ J? exp [— A\J — In E\ , with A > a model 
dependent constant 



50|, |52| . By a real space RG in the strong disorder regime 



66, 



67, 



6jj and 



through a supersymmetric [6J| and replicate d |7Qj random gauge field theory approach as well 
as within the standard non-linear cr-model |7_l[, it has been recently argued that the correct 
asym pto tic expression of the DOS is instead of the form p(E) ~ ~ exp [— A(— In i?) 2 / 3 ] . In 
Ref. |7jj the origin of the disagreement is identified into the existence of an infinite chain 
of relevant operators which are related to moments of the Q-field and which are coupled 
together in the RG equations. 



IX. THE ACTION WITH VECTOR POTENTIALS 

The quasiparticle charge modes, as well as the spin modes when magnetic impurities or 
a magnetic field are present, are not described by the non-linear cr-model (190)1 . which only 
represents the truly massless diffusion modes. Nevertheless, charge and spin conductivities, 
cr c and a s , respectively, can be still evaluated through the stiffness of the corresponding 
modes, although they acquire a mass term. Alternatively, a c and cr s can be determined by 
second derivatives of t 
or to the spin current 

As explained in Appendix B the source field which couples to the charge quasiparticle 
current is the vector potential 



re action with respect to a source field which couples to the charge 

3- 



A c = A s (A T3 S o + AV 3S i), (103) 
where X s is a symmetric matrix in replica space, or, alternatively, 

A c = \ a (A°T 3 a z s + A 1 r 3 (7^ 1 ) ) (104) 
with A a an antisymmetric matrix. On the other hand the spin vector potential which couples 
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to the spin current is given by 



A s = \ s (A°T a z s + A x r Q a zSl ) 



(105) 



or alternatively 

4 = A tt (iVo + iV3 Sl ). (106) 

In the hydro dynamic limit the action in the presence of a vector potential acquires a new 
term which, up to second order is A, is (see Appendix C) 

" a, ..I. v I \ Q+i 

c' v v c L 
for a charge vector potential, where a c is the bare charge conductivity 

S 2 ^f (V fc e) 2 



Sa = — cr c Tr 
32 



VQ + i e -[Q,A c ]) (vQt -i- c [A c ,Q^ - (VQVQ*)] , (107) 



a. 



7TV 



(El + E 



2\2 



(108) 



while for the spin case 



C ™»m 



vg + ^[Q,A s ]j hQi-l[A„Qi]) -(vgvgt) 



(109) 



where a s is the bare spin conductivity, given by Eq.()83|). In the presence of these terms in 
the action, the generating function Z(A) depends now on A 



Z{A) = J DQe~ s ~ SA . 



(110) 



The Kubo formula both for charge and for spin conductivities is recovered by 

' d 2 lnZ d 2 lnZ* 



dA ' 



OA 1 



111) 



A=0 



We can now calculate, through the one-loop expansion, the corrections to the charge and the 
spin conductivities as responses to the source fields. To this end we expand to the second 
order in A the generating function 

1 



DQe 



1 - S 



Z(A) = J DQe~ s ~ SA z 
where Sai and Sa2 are 

s A i = i jJ dR Tr (^Q(R) [Q j (R),A]) 



[.. — Sa2 + -{Sai, 



Sa2 



p 



dRTr([A,Q(R)] [Q*(R),A]) 



(112) 
(113) 
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and t = f = | for the spin case and t = — , f = - for the charge case. 

7TCT 2 ^ 7TCr c ' •> C ° 

Taking, for charge conductivity, the gauge (jl03|) . we can calculate the second derivatives 
of the generating function, 9 and 9 , at one- loop level by expanding Q in terms of 
W and averaging with the quantum weight e~ s . The gauge (jl(J4j) gives the same results. For 
spin conductivity we take the expression ()l()5j) or alternatively (jlOfij) as the vector potential. 
Through Eq. we find the one-loop quantum interference corrections for charge and for 

spin conductivity, which are summarized in the following Table |n] 

TABLE II: One loop corrections to to the spin and charge conductivity in n = replica limit. 





Coset space 


5a s /cj s 


5<7 C / a c 


Yes chiral, Yes T 


U(4n)xU(4n)/U(4n) 








Yes chiral, No f 


U(4n)/0(4n) 





—2g In s 


Yes chiral, magnetic field 


0(4n)/0(2n)xO(2n) 








Yes chiral, spin flip 


U(2n)/U(n)xU(n) 








No chiral, Yes T 


Sp(2n)xSp(2n)/Sp(2n) 


— 2g In s 


— 2g Ins 


No chiral, No T 


Sp(2n)/U(2n) 


—gins 


—g In s 


No chiral, magnetic field 


U(2n)/U(n)xU(n) 








No chiral, spin flip 


0(2n)/U(n) 


glns/2 


^lns/2 



By this procedure we find that the one-loop corrections 5a c /a c and 5<j s /<j s coincide with 
8a /a in the absence of sublattice symmetry. When sublattice symmetry holds, quasiparticle 
charge conductivity may behave differently from spin conductivity, as it happens when time 
reversal symmetry is broken. Notice that, according to Table |HJ quantum interference 
corrections in the diffusive modes influence also the stiffness of modes which are, on the 
contrary, not diffusive. 

X. THE RESIDUAL QUASIPARTICLE INTERACTION 

So far we have dealt with disorder in d-wave superconductors modeled by a BCS Hamil- 
tonian for free Landau-Bogoliubov quasiparticles. However strong correlation is a crucial 
ingredient of the cuprates. Therefore it is important to understand the effects of the resid- 
ual quasiparticle interactions even within the superconducting phase. Besides producing 
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dephasing scattering processes 



73 



the interactions can renormalize both the density of 



states and the stiffness of the spin fluctuations. In this Section we extend our ana! 



53 



ysi s to 



75 



include the quasiparticle interactions, following the original work by Finkel'stein 
for disordered metals. Moreover, we will extend the Finkel'stein model by including the nest- 
ing property, which requires to add interaction amplitudes with (tt, it) momentum transfer. 
Let us first consider the following interaction contributions to the free energy 



\k\<Ck F 



(114) 
(115) 



T Y ^(Pl) 4n.{P^) ia,(pi + k) C^_^(p 2 ~ k), 

|fc|<fc F 

with a and j3 spin indices and n, m and u Matsubara indices while pi, P2 k are the momenta 
involved. T is the temperature. The sum is performed both over momenta and frequencies. 
Indeed since interactions intermix the energies of the particles we should give up to a fixed 
energy description, used in the non-interacting case, introducing a discrete set of Matsubara 
frequencies. 

p+k, <b+v q+k, co+e 



p.v 



q. e 



FIG. 1: Diagram of interaction in particle-hole channel 



These interactions can be rewritten distinguishing the singlet from the triplet channel in 
this way 

r 

T ^2 -^c n {p l )a c n+U} (p l + k)c m {p 2 )(yoC m ^{p2-k) 1 (116) 
|fc|<fc F 

r 

- T ) — c n {pi) ac n+u) {pi + k) c m (p 2 ) vcm-uipi ~ k), (117) 
with T t = r 2 /2 and T s = Ti — T 2 /2. By gaussian integration and using (jSJ) we have 

e T T. §Cn(Pl) °" Cn+Lj(pl+k)c m (p 2 ) a c m _ CJ (p 2 -fc) _ (118) 
dXe -\Y,^{Xo{u)X Q {-u)-X- i {u)X 3 {~Lu))+2iY J V=TT(x^ 
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being X (— to) = Xq(u), X 3 (— uj) = — X 3 (u>) auxiliary Hubbard- Stratonovich fields and 
T = —T s for the singlet particle-hole channel with a = <jq or V = T t with a = a for the 
triplet particle-hole channel. In Eq. (|118|) the positive-negative energy index of the Nambu 
spinors have been extended to label the positive and negative Matsubara frequencies. By 
integrating the fermionic degrees of freedom, the full action including interaction is 

^Tr In (iUetf - UH^U ] + iQ sp + 2ty^jfX Q Ur a t U^ (119) 

where e nm = e n 5 nm = mrTSik with n, m odd integers. U is the unitary transformation, previ- 
ously called T in the non-interacting case. Then, by expanding in terms of the saddle point 
Green's functions, we find new terms in the action that represents the residual interactions 
in the p-h channels, namely 

) + X 3 (u})Tr(r 3 aQ nin+UJ ) . 

(120) 

Finally, by integrating over the auxiliary fields X and X 3 we get the following contributions 
to the free energy for the singlet channel 

ix 2 v 2 

~ T Y1 ~s~ Ts Yl PMQfi.nW 7 * a o) Tr(Q m+UJim Ti (T )) , (121) 

1=0,3 

and for the triplet channel 

T E^ r * E (^(Qn^^TriQ^^na)) . (122) 

Z=0,3 

In the replica space the Q matrices contained in Eq. (|12ip and in Eq. (|122|) are diagonal and 
have the same indices since residual interactions is present at fixed disorder. For convenience 
we will put upper latin indices, like Q ab , to denote replicas. In d-wave superconductors, from 
[T, r 2 Si] = we have 

t"2Si<3t 2 si = -Q (123) 

together with the condition 

C*Q*C = Q, (124) 

where Sj, i = 0,1, 2, 3, are Pauli matricies acting on the signs of Matsubara frequencies. For 
the singlet and To and r 3 components, this means 

QsO,nm = QsO,mn = ~ Q SO, —n—mi 

(125) 

Qs3,nm = ~Qs3,mn = Qs\-n-mi (126) 
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having denned 

Q = Qsvo + iQr-° (127) 

in spin space and 

Qs = QsoTo + i QsjTj, Qr = QToTo + i 2j QrjTj (128) 

3=1,2,3 3=1,2,3 

in particle-hole space. Because of Eqs. (J1251 1126J1 the interaction in the p-h singlet channel 
is therefore 



+ Q S3,n,n+wQ 



aa 

S3,m+w,m 



nmus a 



\ ^ \ r^aa r^aa , r^aa y^j 



aa 

S3,— n— uj,— S3,m+uj,m' 



By setting — n —>■ n + u in the last term, we recover the first with opposite sign, hence the 
sum is zero. More specifically, if we consider the transformation Q n , m — > Q_ m _ n we find 
that the interaction terms only involve the symmetric components of the Q's because of 
energy conservation. On the other hand, only the antisymmetric Qs and the symmetric Qt 
stay massless in the superconducting phase due to Eq. (|123|) . This means that the singlet 
term, with T s , is suppress in ci-wave superconductors by symmetry and only the triplet, with 
T t , survives. This is physically expected being charge fluctuations not diffusive. 

We now take into account also the interaction in the Cooper channel. The diffusive 
cooperon represents fluctuation in the particle-particle channel with s-wave symmetry. Since 
the real part of the order parameter is already finite, fluctuations in the r 2 si channels are 
not diffusive, while only fluctuations in the t%s\ channel, corresponding to fluctuations of an 
is order parameter, stay massless. In the presence of residual interaction in the p-p channel, 
we must also consider the term 

T E ^<{Pi)ct-n{k- Pl )ci{p 2 )cZ_ m {k-p 2 ). (129) 

|fc|<fcF 

By introducing one more auxiliary field, Y a ^ = {Y^ a )* with a and /3 spin indices, the 
quantum weight due to p-p interaction can be rewritten as 

' dYe' 1 ^ y ™"( fc ) y ™ a (- fe )+ i £ v / m(5S(Pi)X?' 3 (p 1 +P 2 )5S_„(p 2 )+cg(pi)y^ a (pi+P2)cs_ n (p 2 )) 

From Eq. © the following equalities hold 

*(n + iT 2 )V = -ic<T y c, (131) 
- w- 2 )W = -ica y c, (132) 
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FIG. 2: Diagram of interaction in particle-particle channel 



so that, calling Yf = £ 7 Y" 7 ^ an d Y f = E 7 < 7 ^ 7/3 , which im P lies Y t = ( Y f)*, E Q- 
(|13Uj) becomes 

J dy e -3 E YgY£+Z VTT- c (KY^r+^- n +^Y^r-^Z- n ) ? ( 133 ) 

where r — T\ ± 2x2. Integrating over fermions we find 

y ^y e -fE ¥v^(^KQ^„,„- + )+^Tr(Q^ mim r-)) ^ 

and finally, after integration over Yr, we obtain the following additional term to the free 
energy, representing the interaction in the Cooper channel 

-^ T cTr spin {Tr(Q n+0J _ n T + )Tr(Q m+LU ,- m T~)}. (135) 

Also in this case the Q matrices are diagonal in replica space and both of them have the 
same replica index. By the charge conjugacy relation C t Q t C = Q, the triplet terms don't 
contribute since 

QTl,nm = ~QTl,mni QT2,nm = ~QT2,mn (136) 

so, if in Eq. (j!35|) we transpose Q n +uj,-n and put — n — > n + to we will have triplet terms 
with opposite sign. This means that at the end only the following term remains 

T E x Fc ^ (Tr(Q:+ ^ - nT; " o) Tr - mTi ffo)) • (137) 

1=1,2 

This result is valid also in metal phase and reflects the local character of the Q matrix. 



XI. RENORMALIZATION GROUP WITH INTERACTIONS 

Now we calculate the corrections to the conductivity and to the density of states due to 
the interactions. Let us first consider the model without sublattice symmetry, for which the 
relevant interactions are T t and r c defined above. 
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A. Interactions with small momentum transfer 



The properties of massless modes in the Matsubara frequency space are the following, 
having imposed antihermitianicity and conditions (|31j) . (|33j) . 
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where a and h are replica indices, while n and m are odd integer Matsubara indices with 
opposite sign. If we take n = —m we recover the symmetry properties derived in the 
previous chapters and reported in Appendix A. By these relations we can write the gaussian 
propagator 



(w£ nm (k)W£ rq (k)) = (±)-{l-X n X m )D nrn (k) 



rac sbd r i 1 rad rbc / \i rac r 

nr mq [ J ^nq ^v^^ ' ' — *■ " 



bd 



i r_ii cad rfec 

mq L~^J "ng "mr V / "n— r ^m—q V / [ J 'n—q ^m—r 1 1 



138) 



where S means S for singlet and T for triplet components, 8™ = 8 ac 8 nT , (±) are related to 
real or imaginary matrix elements of W, [±] for symmetric or antisymmetric matrix, (— ) % 
the sign that W acquires changing the signs of Matsubara frequencies and this occurs only 
for modes proportional to T\ and r 3 , and finally 

1 1 



D nm (k) 



(139) 



Atiu Dk 2 + z\e n - e m \ ' 

with sign(e n ) = X n = —X m where e n = nnT is a fermionic Matsubara frequency with odd 
integer n. The factor z is the frequency renormalization and D = oj{2v) the diffusion 
coefficient. 

Let us introduce slow and fast modes in the spirit of Wilson Polyakov procedure, as we 
have done before, 

Q = U\Q f U s = Utu}Q sp U f U s , (140) 
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FIG. 3: Diagram of the diffusion propagator 



with U = e 2 (without chirality, U = U) and finally Q sPnm = A n 5 nm . 
U s contains only slow momentum fluctuations and 

Usnm = $nm, if (s e r) _1 < |e n | < r -1 or (s e r) _1 < \e rn \ < r^ 1 , (141) 

where r -1 acts as an energy cutoff and s e > 1 is the rescaling factor. The massless fast 
modes satisfy by definition 

W fn Jk) = if {Dk 2 , \e n \, \e m \} < (s.r)- 1 . (142) 

Now let us expand the interaction contributions to the full action, namely Eqs. (|121I122I137|) 
multiplied by T _1 , in terms of Wf, leaving slow U s unexpanded. In this way, besides the 
terms (f§H|) and (JHSj), also the following contributions should be evaluated in the one- loop 
expansion 

2 

S(ni =F n 2 ± n 3 - n 4 ), (143) 

Tfi 2 ±n 3 -n 4 ), (144) 

where the upper indices are in the replica space, T = T t , a = a, and we sum over t\ = r , r 3 
and energies with constraint 5(rii — n 2 + n 3 — n 4 ) for p-h triplet channel and V = T c , a = a , 
ti = ti,t 2 , with energy conservation law 8{n\ + n 2 — n 3 — n 4 ) for p-p Cooper channel (in 
metal phase we would have also the singlet channel with T = —T s , o = <jq, t\ = tq,t 3 and 
energy constraint 5(jii — n 2 + n 3 — n 4 )) and finally Q = WQ sp U is the slow mode matrix. 
The sums in Eqs. (I143[ I144|) are performed over all upper (replica) and lower (frequencies) 
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indices and over /, according to the channel. Before starting the renormalization of all the 
parameters involved in the theory we can easily take into account all the ladders diagrams 
sketched in FigEJ substituting the bare particle-hole triplet scattering amplitude with 

Dq 2 + z\uj\ 



T t (q,uj) = IV 

Dq 1 + [ z + 2pT t )\uj\ 

which is the algebraic infinite ladder summation 77 1. 



(145) 



FIG. 4: Ladder summation 



For Cooper amplitude we have instead, already after summing the first two diagrams, a 
logarithmic divergent term which does not depend on g, the small parameter of the theory. 
To calculate one-loop corrections to conductivity due to the interactions, we have to consider 
the following averages over fast modes 

(3D - (3iSD - (S*SD + \({S 2 fSD, (146) 

where S l and S 2 are defined by Eq. (jHEJ) and Eq. (jHHJ). For instance, contracting the fast 
modes, the contribution coming from 

2 

( S lnt) = -~Y $^r(p,e ni ± e n2 )D mim2 (p + pi)(^l - A mi A m2 ) <J(m =F n 2 ± n 3 - n 4 ) 
[Tr{\ m2 U^ n2 { Vl -P2)r^^ i (p 2 )A mi C/^ n4 (-p 3 )r^^* 2 (p3 

- (-) ; Tr(A m2 ^f 2n2 ( Pl -pa)^* >)A mi <_^(-p 3 )^*^(p 3 -piX)}T) 

which renormalizes the spin conductivity is obtained summing over fast momenta m\ (the 
sum over slow momenta, instead, renormalizes the scattering amplitudes), giving simply 

-n 2 un a Y,^nP,e mi ±e n2 )D mim2 (p + Pl )Tr(ul% 2 ^ (148) 

with n a = 3 in the triplet p-h channel and n a = 1 in the Cooper one. Expanding the 
propagator D in terms of slow pi, neglecting slow frequencies both in T and in D, we have 
a factorized term 

- n a jJ2piTr(u}X(-PiW 9 ^M)) = n ° J j drTr(AA), (149) 
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where A = VUU* 1 (notice that here, since chirality is not present, U = U) and 
J = ±J dpdeT(p,e)^ {Dp2 D + p zlel)2 - ^Sr] 5 )- The other mean values in Ec l- dffl) give 
contributions of the same kind and also proportional to J drTr(AXAX) (with chirality 
J drTr(AX'-fiAX'-fi)) . Reminding that 

2Tr(AXAX - AA) = Tr(VQVQ ] ), (150) 

we obtain corrections to the stiffness coefficient. 

To calculate corrections to the density of states, besides (S u ), we will consider 

(SLSu), (151) 

where 

S » = \ E Tr (^Il KmW^W^UfcJ. (152) 

In order to calculate corrections to the interaction amplitudes at first order in g, we need 
to consider the sum 

(Sint) — -({Sint) 2 ) + ^((Sint) 2 Si nt ) + -^(Si nt (Sf nt ) 2 ) — — ((Sl nt Sf nt ) 2 ) , (153) 

where S int = S} nt + Sf nt . Because of the different structures of the two terms Sj nt and Sf nt , 
we notice that (Sf nt ) in the first term, (S} nt Sf nt ) in the second term and the third term of 
Eq. ()153j) are vertex contributions (in FigEH the corresponding diagrams are sketched using 
only the p-h triplet scattering amplitude) while the last two in Eq. (|153|) and ((Sf nt ) 2 ) in 
the second term are bubble corrections, (see Fig. |7J). In the normal phase the the first two 
diagrams of Fig. [7| and the first diagram of Fig. |U1 vanish. There are other diagrams which 
give corrections to T t in which r c appears and viceversa, with the same topology but with 
different energy-momentum conservation laws. 

Finally, the renormalization of z is determined by three contributions, the first two result 
from the corrections to the density of states (S u ) and from Eq. ()151j) and the third from the 
expansion of {S} nt ), written explicitly in Eq. (|148J) . in terms of the slow frequency n 2 . This 
latter contribution cancels exactly all the terms of order higher than one in the interaction 
amplitudes coming from Eq. (|15ip . 
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FIG. 6: Vertex contributions. 




FIG. 7: Bubble contributions. 



B. RG equations without chiral symmetry 

The final renormalization group equations for ci-wave superconductors, including residual 
quasiparticle interactions T t and T c , in all the cases collected in Table are the following: 
a. No chiral, Yes T. If time-reversal symmetry is preserved we obtain 

dg/dl = g 2 {-3uT t f 2 {z, z, z 2 ) - vY c jz + 1} , (154) 

dz/dl = g{3uT t + uY c - z/2} , (155) 

duT t /dl = g {v (r t + T c ) /2 + 4z/ 2 (r t r t + T t T c ) /z + 4z/ 3 (T t T c T t ) /z 2 - uT t - u 2 T t T t ^) 

dvTjdl = g{v (3r t - T c ) /2 + 2z/ 2 (3r t r c *A(*, z 2 ) + T c r c ) /z - vT c ) - 2(vT c ) 2 /z, (157) 
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where, as usual, we have introduced the functions f\ and f 2 which result from integrations 
of products of diffusion propagators, defined by 



A (a, b) = m(a/6)/(a - b), (158) 
f 2 (a,b,c) =2{bf 1 (a,b)-cf 1 (a,c))/(b-c), (159) 

and z 2 = z + 2^r t . Notice that in Eq. ()154j) the last term is taken from Table H] with n = 
and divided by 2 since here we have adopted an energy scaling factor instead of a momentum 
scaling factor, being I — lns e with s e defined in Eq. ()141j) . The contributions due to the 
non- vanishing density of states in the non-interacting case are singled out, they are the last 
two terms of Eq. (J156j) (corresponding to the first diagram of Fig. |U] and to the first two of 
Fig. |7J) and the last term in the curly brackets of Eq. (|157|) . In the normal phase, they 







53. 751. The last term 



disappear meanwhile the singlet particle-hole channel is turned on 
of Eq. ()157|) is due to the ladder summation. The scaling behavior of the density of states 
becomes 

du/dl = gv {3uT t f l (z, z 2 ) + vTjz - 1/2} . (160) 

Notice that the theory has only four parameters, g, z, vTt, vT c since the density of states 
appears always multiplied with the scattering amplitudes. Summing and writing in terms 
of u t = 2i/T t /z and u c = 1vY c jz we can separate from the set of equations the ^-equation, 
obtaining 



dg/dl 


= g 2 {3(1 - (1 + u t ) ln(l + u t )/u t ) - u c /2 + 1} , 


(161) 


dz/dl 


= gz{3u t + u c - l}/2, 


(162) 


dut 1 dl 


= g {u c /2 + 3 u t u c /2 + u t u c u t } , 


(163) 


du c / dl 


= g {(3u t - 2u c ) /2 + 3 (u c ln(l + u t ) - u t u c ) /2} - (u c ) 2 , 


(164) 


du/dl 


= gu{3ln{l + u t ) + u c -l}/2. 


(165) 



Supposing T -1 , the inverse of the temperature, being the coherence time for the quasipar- 
ticles, the integration of the equations above will run from T to r _1 , the elastic scattering 
rate. Strictly speaking, half of the localization correction in the last term of Eq. ()161|) . is 
cut off not by the temperature but by the inelastic scattering rate r 4 ~ . However, r^ 1 , for 
some source of inelastic processes, it is found to be linear in T . At the first order in the 
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scattering amplitudes, supposing them so small, at the moment, that they do not flow fast, 
we have the following variation in the conductivity induced both by the disorder and by the 
interactions 



5a 



1 



1 



u c ln(7Y) 



(166) 



2tt 2 V" 2"" 2 

Solving now the whole set of equations (jl61H164j) we find that, starting with a small 

and positive u t o (in FigJE] u t o = u t (T ~ r _1 ) = 0.15, u c o = u c (T ~ r _1 ) = and 
go = g(T ~ r- 1 ) 



0.03), g increases for the presence of the quantum interferences but 



when the interactions become strong enough, after reaching a maximum, it decreases at low 
temperature, as shown in Fig. |HJ The temperature of crossover depends on the interaction 
amplitudes. The smaller is the value of the positive Uto, the higher is the peak of g, while Ut 



diverges as in the normal phase 
the localization. 



75 
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A negative starting value of w i0 , instead, reinforces 
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FIG. 8: Scaling behavior of g of paragraph a, in the inserts ut and u c . 



b. No chiral, No T . If time reversal symmetry is broken, r c disappears from Eqs. (jl541 
115611160]) and the last term of Eq. ()154|) . which is the contribution due only to the disorder, 
is halved (see Table QJ). u t remains unrenormalized at first order in g fEq. (J163l) with u c = 0), 
namely du t /dl = 0(g 2 ), which is a result obtained also by Keldysh technique 78], so u t = u t o 
at the one-loop level. The equation for g, exact in the interaction, becomes simply 

dg/dl = g 2 {3 (1 - (1 + hi(l + u m )/u m ) + 1/2} , (167) 

implying that there is a positive critical value of Uto {u^ ~ 0.37) for which dg/dl = at 
lowest order in g. Above u^q the resistivity g decreases, below wJq it increases. 
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c. No chiral, magnetic field. If also the spin rotational invariance is broken, we have 



dg/dl = -g 2 vT t f 2 (z,z,z 2 ), (168) 

dz/dl = guT t , (169) 

duFt/dl = -guTt/2, (170) 

du/dl = guT t f 1 (z,z 2 ). (171) 

Also in this case, as in the previous ones, we can eliminate z introducing u t , obtaining 

dg/dl = g 2 {\-(\ + ut) ln(l + u t )/u t } , (172) 

du t /dl = -gu t (l + u t )/2. (173) 



The solution is a transient since the scattering amplitude flows to zero and the conductivity 
changes value only in the meanwhile. For a small value of u t0 = u t (T ~ t _1 ), we can retain 
only lowest order in u t so that, given g = g(T ~ t _1 ), the renormalized resistivity is given 
by 

9 = 9oe~ Uto . (174) 

d. No chiral, magnetic impurities. In the presence of magnetic impurities there are 
no corrections due to the interactions. Indeed both cooperons and triplet diffusons become 
massive. 

C. Interactions with (tt,tt) momentum transfer 

At the nesting point where the staggered fluctuations become diffusive (i.e when the 
chiral symmetry holds), also quasiparticle interactions with (n,n) momentum transfer must 
be included in the effective action. First notice that in the sublattice representation, the the 
contribution due to the interactions with small momentum transfer can be rewritten as 

2 

E ^ vT ° Tr (QZn 2 n^o)Tr(QZ n4 n^ ) S( ni Tn 2 ±n 3 - n 4 ), (175) 

where 70 is the identity in the sublattice space and r° can be — F°, F° t and T° c . 

Let us now consider interactions whose transferred momentum is = (±7r, ±7r), which also 

involve quasiparticles at the Fermi energy in the case of half filling, 

r 3 

T E y^(p 1 )c^(p 2 )c^_ a; (p 2 -A;- g7r )< + > 1 + A; + g7r ), (176) 
|fc|<fc F 
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r 

T ~~2 ^(P 1 ) CmiPz) Cn+wiPi + k + q w ) c^_ w {p2 -k-q n ), (177) 

|fc|<fcF 

r 3 

T Yl ^<^)t-n{k-p 1 + q^ci{ P2 )C_ m {k-p 2 + q,). (178) 
|fe|<fc F 

3 

In this case T indicates the amplitude related to staggered modes. Repeating the steps 
which led to Eq. (fTZT]). Eq. (|T22l and Eq. (fT37)l . we obtain the following additional 
contribution to the action 

2 

£ rr(g- ri2 r^ 73 )Tr(g- n4 r^73) <y(m =F n 2 ± n 3 - n 4 ), (179) 

3 3 3 3 3 3 

with T taking the values — T s = T 2 /2 — I\ for the singlet p-h staggered channel, T t = T 2 /2 

3 

for the triplet and T c for the p-p Cooper staggered channel. In the superconducting phase 
only T 3 and r 3 are relevant. The properties of W 3 massless modes in energy space derive 
from the conditions (JH2J), p3|) and antihermitianicity, leading to 
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where now n and m have same signs. By these properties we can write down the propagators, 
reported in Appendix D, and evaluate Eq. f)146|) . Eq. (j!51|) and Eq. ()153j) . 

D. RG equations with chiral symmetry 

The renormalization group equations in different symmetry classes in the presence of 
chirality in the sublattice space are the following: 
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e. Yes chiral, Yes T. If both time reversal symmetry and spin rotation invariance 
are preserved, we have, together with T = g/c and dc/dl = —4c 2 (see Eqs. (jlOOl I10l]0 . the 
following equations 

dg/dl = g 2 {-3vT° t f 2 (z , z , z 2 ) + 3vY\jz - vV°Jz + uT 3 J 2 (z, z,z c )^ , (180) 
dz/dl = g |3i/r° - 3uT 3 t + uT° c - vV] + zT/8~j , (181) 

duvl/di = g {u (r; + r; + r: + rj) /2 + 4, 2 (rjrj + - r t 3 r t 3 - r° t r 3 t ) ,z 

+ 4,3( r ;r;r;-r;r;r;)/^ 

+rvr t °/4 + rv 2 r^/4* + ivr, 3 (1 + 2vv]/z) 2 /A , (182) 

dvv]/di = ^{z,(r: + r: + r° c + r:)/2-2z, 2 (r° t r 3 t zh(z,z 2 ) + 3r 3 t r 3 t -r° c r 3 t 

+vy t zh(z, Iz + rv(r t 3 + r t °)/4} + 2(uv] f/z } (iss) 
duvl/di = g [v (sr t ° + sr; - r: - r 3 c ) /2 + 2^ 2 (ar^/^, z 2 ) - 3t]t° c 

+Vy c - V 3 f c zh{z, z c )) Iz + I>(l1 + r 3 c )/4} - 2(<) 2 /z, (184) 

dur 3 c /di = g {v (3^ + 3^ - r° c - r 3 c ) /2 + 4z, 2 (r^ - srK) /z+ 



4, 3 (sivivr; - r 3 yy) /z 2 

+TuT :i c /A - Tu'TyjAz + TuT° c (l - 2</,) 2 /4I 



(185) 



dz//rf/ = gv{3vr t f 1 {z,z 2 )-3vr 3 t /z + vr° c /z-vr 3 J 1 {z,z c ) + zr/8y (186) 

3 

where z 2 = z + 2i/T f and z c = 2 — 2z/T c , which come from ladder summations in the triplet 
p-h channel, Eq. ()145j) . and in the staggered p-p Cooper channel, 

r c (g^) = r c ^ 2 + ( ^ 2<)M - (187) 

The last term of Eq. (jl83|) comes also from the ladder summation in the staggered p-h triplet 
channel that present the same log- divergence of the standard Cooper channel, Eq. (jl84|) . 
namely it does not depend on g. As a result, for repulsive interaction, the Stoner insta- 
bility towards a spin-density wave is not destroyed by disorder. At the first order in the 
scattering amplitudes, introducing u° = 2T t /z, u° c — 2Tjz, u? = 2T t /z and u\ = 2Tjz the 
conductivity and the density of states are corrected in this way 

8a /a = g [3« - u 3 t )/2 + (u° - u 3 )/2] m(Tr)-\ (188) 
8v/v = g [3K> - ul)/2 + « - ul)/2 + r/8] ln(Tr)- 1 . (189) 
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f. Yes chiral, No T. If time reversal invariance symmetry is broken both T c and T c 
disappear in the equations above, dc/dl = —2c 2 , and 

dg/dl = g 2 v {-3T° t f 2 (z,z,z 2 ) + 3T 3 t /z} /2, (190) 
dz/dl = gv^3Y\-3T 3 t + z(Y/A-l)/2^ , (191) 

dur t /di = g {u (r t ° + rj) /2 + Av 2 (vy t - t\t\ - rjrj) /z -aStXtI/z 2 

+(r/4 - i)vv] + (r/4 - i)v 2 vy t /z + YvY 3 t (1 + 2uy]/z) 2 /a} , (192) 
dvV]/dl = g{u(r° t +^) /2-2v 2 {T t r ) t zf 1 (z,z 2 ) + 3T s t r i ) /z 

+ (r/4 - l)vV] + TuT° t /A} + 2(uT 3 t y/z, (193) 
dv/dl = gu^3uT° t f 1 (z,z 2 ) -3vY]/z+ (r/4- l)/2 j> . (194) 

g. Yes chiral, magnetic field. If the c?-wave superconductor is embedded in a constant 
magnetic field, namely the spin rotation invariance is broken, we have 



dg/dl 


= g 2 {-vF° t f 2 (z, z, z 2 ) + 2vY]/z - 1/2 } , 


(195) 


dz/dl 


= g{vT° t -2vr 3 t }, 


(196) 


dvfjdl 


= g (-rJ/2 + rj) - 4^ r ; r ;/,} , 


(197) 


dvY\/dl 


= g {vV]/2 - 2z/ 2 (3r t r^A(^, z 2 ) + 2TX) / z) + 2(vT]) 2 /z, 


(198) 


dv/dl 


= gv{vT° t f 1 {z,z 2 )-2vT]\. 


(199) 



An effect of the SU(2) symmetry breaking is to reduce the triplet scattering amplitudes as 
r° — > r°/3 and T'l — > 2^/3 in the conductance and in the DOS. Indeed for the standard 
modes the massless ones are those which commute with the Zeeman term, i.e. only one 
component in the triplet sector survives, while the staggered massless modes anticommute 
with the Zeeman term, i.e. there are two components in the triplet sector. Solving the 
equations above we obtain that for moderate values of the interactions g decreases, w° is 
finite and wf does not diverge if the starting value u^ ~ — g 1 ! 2 or if u° is positive and u^ ~ 0. 
h. Yes chiral, magnetic impurities. In the presence of magnetic impurities we have 



dg/dl = g 2 vTJ 2 (z,z,z c ), (200) 

dz/dl = gvT 3 c , (201) 

dvT 3 Jdl = -gvT 3 c /2, (202) 
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dvjdl = guT c f 1 (z,z c ). 



(203) 



These equations are formally the same of Eq. (|168H171|) . consistently to the fact that the soft 
modes live in the same manifold, U(2n)/U(n)xU(n), where now n = ^replicas x ^positive 
frequencies. At low interaction regime we have the solution g = g e n *, similar to Eq. ()174|) 
except for the sign in front of the interaction amplitude. 

E. In normal phase 

The role of staggered fluctuations can also be considered in the metallic phase, extending 
the Finkel'stein model to include nesting effects. If superconductive order parameter is 
turned off and if staggered fluctuations are supposed to be massless, the singlet contributions 
stay relevant and the renormalization group equations become the following: 

i. Yes chiral, Yes T. In the presence of time reversal symmetry we obtain 

dg/dl = g 2 ^uT° a f 2 (z,z,z 1 ) - uT 3 Jz - 3uT° t f 2 (z, z, z 2 ) + 3vV\jz - 2uT° c /z + 2uT 3 J 2 (z, z, z c )^ , 
dz/dl = g |.r 3 s - uT" s + 3^r° - 3vV] + 2uT° c - 2vT\ + z (1 + r/8)} , 

duvl/di = g \u (r; + r; + r t ° + rj + 2r: + 2^) /2 + 4, 2 (rtf + r^ + 2^ - rrrl 

-ry t ) iz + 4,3 ( 2r ;r: r ; - rjrjr, + r ; r ; r ;) iz 2 

+ (2 + r/4) (uT° t + uXK/z) + TuT't (l + 2uT° t /z) 2 , 

duvl/di = ^{,(r: + r: + r; + r; + 2r: + 2r:)/2-2z. 2 (rXzf^z^ + srX 

-2vy t +2vy t zf 1 (z,z c )+vy t zf 1 (z,z 1 )-vy t ) / z 

+ (2 + T/A)X+TuT t /A}+2(url) 2 /z, 

dvvl/di = g {u (-r° s - r; + sr: + sr; + 2r: + 2^) /2 + ^ (?y s - 3r 3 t r° s 
+rX + vXz/z) iz + 4,3 {-vyy s + ar^rj) /z 2 
-4, 2 (2r;r:/, c + r;r;(I-I)) 

+ (2 + r/4) (uv° s - StX/z) + r^ (1 - 2uY°jz) 2 , 

dufjdi = g {u (-r° s - r; + srj + X + 2V° C + 2^) /2 + 2, 2 (-ry s zh(z, Zl ) 
+ry s + 3T°y s zf 1 (z, z 2 ) - 3vy s + 21^ - 2rK*/ 1 (* l z c ) 
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+4r° c r^/ 1 (z, z c )) /z+{2 + V/A) uV] + VvV°Ja) - 2(uT 3 s ) 2 /z, 

duT°jdi = g [u (r: + r; + sr? + sr, 3 ) /2 + 2^ (3r t °r^A(z, z 2 ) - sr 3 ^ + r 3 s r° c 
-r°r°z/i(z,zi) + 2r°r° -2r"r°z/i(«,« c ) + 2r°r°«/i(z,ari) 
+2r:r>/ 1 (^, z c )) /* + (2 + r/4) v r° e + rvrI/4} - 2(<) 2 /z, 

durl/di = g [ v (v° s + r 3 s + ar t ° + srj) /2 + Au 2 (tX + rXzWzc, Zl ) - h(z, Zl )) 

-3TX + r 3 s V° c ) /z + Au 3 (3Vy t Tl - 21^ + 1^) jx 



>z 2 



i 2 r|T; ( A^x, z c ) + Zl ) - 2f,(z, z c ) + 1 



+ (2 + T/4) (uY 3 c - u*rX/z) + TuT° c (l - 2uI fl Jz) * , 

d!//dz = ^z/ 2 jr 3 - r° s zft(z, Zx ) + srlzhiz, z 2 ) - 3r* + 2T° C - 2r^/ 1 (*, z c ) + z (i + r/8)} /, 



where z\ = z — 2uT s and z%, z c as before. Z\ is related to the compressibility 76], here it 
appears in the diffusion propagators, modified by singlet p-h channel ladder summation 

o . o Dq 2 + z\u\ 
1 s W> ^ — 1 s n 2 i / o — f' 

If chiral symmetry is broken, the last contributions in the curly brackets proportional to 
(2 + T/4) and to T together with all the contributions proportions to the T 3 s vanish and we 

n t 

obtain again the standard Finkel'stein equations |53J,175|. The equations above, although 
rather complex, hide a very amazing property, they in fact are symmetric with respect to 
the following transformation 

r° r 3 

r 3 r° r 3 

r , = r c < — > -r t . 

This symmetric property represents the invariance with respect to the particle-hole symme- 
try transformation, = d^, = (— ) l d\^ that maps the charge to the spin and viceversa. 

j. Yes chiral, No T. If time reversal symmetry is broken the renormalization group 
equations are found to be formally the same obtained in the superconductive case, when 
time reversal symmetry holds, Eqs. (jl8()H18fi)) . with the substitutions 



r! - r! 



r! - r° 
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This accidentally simple equivalence of the equations is consistent with the fact that in both 
the cases the soft modes take values in the same coset U(4n)xU(4n)/U(4n). 

k. Yes chiral, magnetic field. In the presence of a constant magnetic field we have 

dg/dl = g 2 {vT° s f 2 (z, z, zj - uT° t f 2 (z, z,z 2 ) + 2vY\/z^ , 

dz/ dl = g | -vT° s + vT° t - 2uV\ j , 

dvY°Jdl = g{u {t°J2 - r?/2 + v]) - torfft/z) , 

dvYl/dl = g{vV\/2-2v 2 (tfftzf^z, z 2 ) + 21^ + Y° a Y\zh{z, z x )) / z) + 2{vV])\ 

dvfjdi = g{u(-r: + rl + 2rl)/2-8u 2 r:r:/z + Au 3 ry t r° s /z 2 }, 

dz/dl = 01/ j-i/r°/i(z,2i) + uT° t fi(z,z 2 ) - 2z/r t 3 | . 

I. Yes chiral, magnetic impurities. While with magnetic impurities we have 

dg/dl = g 2 {uT° a f 2 (z, z, z x ) + 2vT*Mz, z, z c ) + 1/2} , 
dz/dl = -g^uT s +2uT 3 c + z/2\ , 

duTl/dl = g\vTl - \vT a s + u 2 Vy s /z + Av 2 tX/z c - Au 2 (21^/* + T° S T° S (1 - I 
dvTl/dl = g^T^ + ^Ty^hiz^z^-hiz^z^-vTl + ^Tyjz 

-4z/ 2 r'r' U(zi, z c ) + h(z, z x ) - 2A(z, z c ) + ~ 
dv/dl = -gu[uY°J 1 {z,z 1 ) + 2uYlf l {z ) z c ) + 1/2]. 
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Notice that in this case by the transformation r s = T c — * — r t one obtains the same 
equations of the superconductive case with broken time reversal and chiral symmetries, 
analyzed in paragraph b, since also in the latter case the soft modes take values in the same 
manifold Sp(2n)/U(2n). 

Now we can also study the scaling behavior of the quantities z\ = z — 2uT° s and z 2 = 
z + 2vY\ that are related b y W ard identities to the compressibility and to the static spin 
susceptibility 

H- As in 00, because the ladder terms do not involve the small parameter 

3 

g, we can consider the situation in which only the lowest order with respect to T c and T t 
(by choosing T° c0 > g 1 ! 2 and T 3 t0 < —g 1 ^ 2 ) are retained in all the equations, consequently, we 
have for z\ and z 2 the following equations: 
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(z) if time reversal symmetry holds 

d (z - 2uT° s ) /dl = 2gv ^ - 3 rJ - 21^ + Z - [\ + Q - E^) (s - 2^ 

d (* + 2<) /dZ = 2gv fr 3 s - r t 3 + 2T° C + 2T° t ^ + X - (l + E) + E^) (z + 2<) " /, 2 

(j) if time reversal symmetry is broken 

d(z- 2uT°\ /dl = 2gv {v\ - 3r" + zT/lQ - VT 3 JA\ (z - 2vY°\ 2 /z 2 , 
d(z + 2vrty /dl = 2gv (v\ - v\ + 2V\^- + zT jVa + (z + 2uY° t ^ / z 2 



,2 



(fc) with constant magnetic field 



2 

2 



d (z - 2vr°) /dZ = -4guT 3 t (z - 2uT° s ^ /z 
d(z + 2url) /dl = 0, 

(/) with magnetic impurities 

d(z-2uT° s )/dl = -2gu[2Yl^ + -^ (z - 2^ / z\ 

Notice that in all the cases, when chirality is spoiled, z\ is not renormalized as already 
pointed out previously in Refs. 

XII. CONCLUSIONS 

In this work we have analyzed the role of disorder in d-wave superconductors, which have 
gapless Landau-Bogoliubov quasiparticle excitations. We have considered several universal- 
ity classes, including the chiral symmetry which occurs at half-filling for a two-sublattice 
model. In addition, we have studied in details the effects of the residual quasiparticle inter- 
action. The main results of this work are summarized in the following. 

In the presence of non magnetic impurities the spin conductivity is suppressed by quantum 
interference corrections, in agreement with Ref. |32|. The density of states vanishes in the 
insulating regime. 

On the contrary, surprisingly, magnetic impurities gives a derealization correction to the 



conductivity meanwhile enhancing the density of states 
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If chiral symmetry is present, namely at half filling for a two-sublattice model, the spin 
stays delocalized in spite of disorder and the conductivity remains finite. The DOS diverges 
in the absence of magnetic field and magnetic impurities 41 1. 

The quasiparticle charge conductivity, namely the optical conductivity at small frequency, 
has in general the same behavior as the spin conductivity. However, when chiral symmetry 
holds and time reversal symmetry is broken, the dirty <i-wave superconductor behaves like a 
spin metal but charge insulator (at finite small frequency, excluding the Drude peak of the 
superflow), manifesting a sort of spin-charge separation. 

Charge fluctuations as well as fluctuations of the real part of the order parameter, as- 
suming the average value to be real, are not diffusive in a superconductor. Therefore the 
residual quasiparticle interaction written in terms of the diffusive modes only contains the 
spin-triplet particle-hole channel and the Cooper channel representing s-wave fluctuations 
of the imaginary part of the order parameter. 

For repulsive interaction, particles and holes repel each other in the spin-triplet channel, 
hence opposing localization. In fact we find that a repulsive residual interaction gives a 
delocalizing correction to the conductivity and enhances the density of states. Namely the 
interaction can compete with quantum interferences due to the disorder. 

According to the symmetry of the system, we have found that the interactions can be i) 
relevant, when both time reversal symmetry, T, and spin rotational invariance, SU(2), are 
preserved, ii) marginal, when T is broken and SU(2) holds, iii) irrelevant, when both T and 
SU(2) are broken. 

We have also studied (jr, tt) momentum transfer interactions, since they are coupled to 
diffusive staggered spin fluctuations at half-filling. We find that the corrections to the con- 
ductivity due to the interaction at (tt, tv) have opposite sign than the corrections coming from 
the interaction at small momentum. Moreover, we notice that an analog of the Anderson 
theorem for s-wave superconductors holds at half-filling for staggered density fluctuations. 
Namely, since these modes are diffusive, the staggered susceptibility remains log-divergent 
even in the presence of disorder. As a result, for repulsive interaction, the Stoner instability 
towards a spin-density wave is not destroyed by disorder. 

In summary we have derived the complete sets of RG equations for many different sym- 
metries at the lowest order in the small parameter g, which does not depend on the disorder, 
and at all orders in the interaction amplitudes, uT t , uF c , which are the real phenomenolog- 
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ical parameters, related instead to the impurity concentration through the density of states 
proportional to the scattering rate, v oc S = We have derived the RG equations even in 
the normal phase when nesting effects occur. In the latter case we find that the compress- 
ibility is a renormalized quantity unlike the standard case [3, 0] when the chiral symmetry 
is not present. 
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APPENDIX A: TRANSVERSE MODES 

In this Appendix we summarize the properties of the W matrices for soft modes needed 
to evaluate the gaussian averages. S (A) and R (I) mean symmetric (antisymmetric) in 
the replica space and real (imaginary) s,, i = 0, 1,2,3, components defined in Eqs. (|53H5Bj) . 
Notice that the number of soft modes decreases by decreasing the symmetry since some of 
them become massive. 

1. Without sublattice symmetry 

(a) Without superconducting order parameter 

i. With time reversal symmetry: Sp(4n)/Sp(2n) x Sp(2n). 



W 


Si 


S2 


W so 


A, R 


S,I 


Wsi 


A, I 


S,R 


W S 2 


A, I 


S,R 


W S3 


S,R 


A, I 


W T0 


S, R 


A, I 


W Tl 


S,I 


A, R 


W T 2 


S,I 


A, R 


W T 3 


A, R 


S,I 
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ii. Without time reversal symmetry: U(4n)/U(2n) x U(2n). 



iii. With magnetic field: U(2 



Wo 


Sl 


S2 


W so 


A, R 


S,I 


W S3 


S,R 


A, I 


Wto 


S,R 


A, I 


W T 3 


A, R 


S,I 


i) x U(2n)/U(2n) 


W 






Wso 


A, R 


S,I 


W S3 


S,R 


A, I 


W Tz0 


S,R 


A, I 


W Tz 3 


A, R 


S,I 



iv. With magnetic impurities: U(2n)/U(n) x U(n). 



W 




S2 




A, R 


S,I 


W S 3 


S,R 


A, I 



(b) With superconducting order parameter 

i. With time reversal symmetry: Sp(2n) x Sp(2n)/Sp(2n). 



Wo 




S2 


Wso 


A, R 




W S1 




S, R 


W S 2 


A, I 




W S 3 




A, I 


Wto 


S, R 




W T i 




A, R 


W T2 


S,I 




W T 3 




S,I 
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ii. Without time reversal symmetry: Sp(2n)/U(2n). 



Wn 
' v 








A, R 




w S3 




A, I 


— * 


S,R 




W T 3 




S,I 



iii. With magnetic field: U(2n)/U(n) x U(n). 



Wo 


Si 


S2 




A, R 




^S3 




A, I 


W Tz0 


S,R 




W Tz 3 




S,I 



iv. With magnetic impurities: 0(2n)/U(n). 



w 


Si 


s 2 


Wso 


A, R 




W S3 




A, I 



2. With sublattice symmetry 

(a) Without superconducting order parameter 

i. With time reversal symmetry: U(8n)/Sp(4n). 



W 


si 


S2 


W so 


A, R 


S,I 


Wsi 


A, I 


S, R 


W S2 


A, I 


S, R 


W S3 


S, R 


A, I 


W T0 


S, R 


A, I 




S,I 


A, R 


w T2 


S,I 


A, R 


W T 3 


A, R 


S,I 





so 


«3 




S,I 


S,I 


Wsi 


S, R 


S, R 


W S2 


S,R 


S, R 


W S 3 


A, I 


A, I 


W T o 


A, I 


A, I 


W T i 


A, R 


A, R 


W T2 


A, R 


A, R 


W T 3 


S,I 


S,I 
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ii. Without time reversal symmetry: U(4ra) x U(4n)/U(4n). 



W 


Si 


S 2 


Wso 


A, R 


S,I 


W S3 


S,R 


A, I 


w T0 


S,R 


A, I 


w T3 


A, R 


S,I 



w 3 


so 


s 3 




S,I 


S,I 




A, I 


A, I 


Wro 


A, I 


A, I 


^T3 


S,I 


S,I 



iii. With magnetic field: U(4ra)/U(2n) x U(2n). 



W 


Si 


s 2 




A, R 


S,I 


^53 


S,R 


A, I 




S,R 


A, I 


W Tz3 


A, R 


S,I 



iv. With magnetic impurities: Sp(2n)/U(2n). 



W 


Sl 


s 2 


Wso 


A, R 


S,I 


w S3 


S,R 


A, I 



(b) With superconducting order parameter 



Wo 


Sl 


s 2 


Wso 


A, R 




W S1 




S, R 


W S2 


A, I 




W S3 




A, I 


W T o 


S, R 




W Tl 




A, R 


w T2 


S,I 




W T 3 




S,I 



W 3 


s 


s 3 




A, I 


A, I 




S,I 


S,I 


W Ty0 


A, I 


A, I 


W Ty 3 


S,I 


S,I 


l). 


^3 


s 


s 3 




S,R 


S, R 




S,R 


S,R 


U(4n)/U(4n). 


^3 


s 


s 3 


^50 


S,I 








S, R 


^52 


S,R 




^53 




A, I 




A, I 








A, R 




A, R 




W T 3 




S,I 
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ii. Without time reversal symmetry: U(4n)/0(4n). 





Sl 


S 2 




W 3 


so 


S3 




A, R 








S,I 




W S3 




A, I 




W S3 




A, I 




S, R 








A, I 




w T3 




S,I 




w T3 




S,I 



iii. With magnetic field: 0(4n)/0(2ra) x 0(2n). 



W 




S2 




A, R 




^53 




A, I 




S, R 




WV z3 




s,i 



iv. With magnetic impurities: U(2n)/U(n) x U(n). 



w 


Sl 


s 2 


Wso 


A, R 




W S 3 




A, I 



W3 


s 


s 3 


W Tx o 


A, I 




W Tx 3 




S,I 


W Ty o 


A, I 




W Ty3 




S,I 


x U(n). 


^3 


so 


S3 






S, R 


W S2 


S,R 





APPENDIX B: CHARGE AND SPIN VECTOR POTENTIALS 



A 



u 



Let us consider an operator diagonal in the Nambu space. Namely 

A tlJ 
A Uij 

where the matrix elements are matrices in the retarded/advanced and replica space. If we 
take An = 0, then such an operator corresponds to 



In the path integral formalism, a generic operator diagonal in the Nambu space, 
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with 




corresponds instead to 

ij 

By comparison we have that 

[Alji + Vy^ijVy] = [(A htJ - A\ :ji ) + a z (A Uj + A f Ui )] . 

Suppose that the operators in question are currents. Then Aij = —Aji, and the above 
relation reads 

[-Alij + o y A 2 ^a y \ = [(Aw + A\ ii3 ) + o z (A Ui - A\^)] . 

In general we can consider either a charge current, implying A] = A[ = A, or a spin current, 
in which case A-\ = —A± = A. 
In the former case 

+ °vMrt<r y ] = [{Aj + A%) + a z (Aij - Alj)] , 

while in the latter 

[-Alij + <? y A 2 , l3 <r y ] = [(Aj - Alj) + a z (Aj + A%) ] . 

We therefore see that, if A (we assume the same property holds for A) is a symmetric 
matrix, the charge current operator is proportional to the identity in spin space and —Ai ti j + 
A 2 ,ij = 2 A^, namely A 2ji j = —A lti j = Aij, while the spin is proportional to a z and —A lti j + 
cr y A 2:ij a y = 2a z A^, implying A 2jij = A hij = -a z Aij. 

In the opposite case of an antisymmetric A, the charge current multiplies o z and 
[Ai ti j + <j y A 2) ij<jy\ = 2a z Aij, leading to A\^j = —A 2ti j = o z Aij, while the spin is propor- 
tional to the identity and [A\^j + o y A 2 ^j<j y \ = 2 Aij, leading to A± ti j = A 2>i j = Aij. 

These relations imply that the charge current , if A is symmetric, is the identity in spin 
space, otherwise is proportional to <r z . For the spin current, the opposite occurs. Let us now 
see which are the vector potentials for charge and spin. 
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1. Current-current correlation function 

Let us suppose to calculate the current- current correlation function (J(R)J(R')}. The 
current in Nambu spinor representation is 

J(R) = -iJ2(R-Ri)^R^RRi^Ri = *Ri J rM R )***> (Bl) 
Ri R1R2 

with 

— * — * 

Jr 1 r 2 (R) — —i{Ri — R2)'Hr 1 r 2 8r 1 r. (B2) 

In the case of spin current we have 7i = (t + iAr 2 s 1 )a z , while for the charge current H = tr 3 . 
The correlation function becomes 

(J(R)J(R')) = (^r, j rMR)^r^RsJrsR,(R')^iu)- (B3) 

R1R2R3R4 

Let us introduce for simplicity multilabels indicating replica, hole-particle, spin, energy and 
position indices. The correlation function can be rewritten in this way 

(¥j^¥j lm m m ) = -J %J (^¥)J lm (m m ¥) + J tJ (WV m }J lm (¥¥}, (B4) 

or in matrix language in the following way 

(fjMjtf) = -Tr (J(M) J(M)) + Tr ( J(^¥) J* . (B5) 

Since 

* = (Ctt)*, (B6) 
with C = iriCTy, the correlation function becomes 

- Tr ( J(M) J(M) + J(M) CJ*C"(##)) . (B7) 
In terms of single particle Green's functions 

- Tr (JG{J + CJ'C^G) . (B8) 
Let us define the following quantities, for spin 

jS Rl R2 = -KRi - R2){t RlR2 + tA RlR2 r 2Sl ), (B9) 

and for charge 

jR lR2 = -KRi-R2)t Rl R 2 , (BIO) 
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such that J = J s a z for spin and J = J c t 3 for charge. In both the cases the following 
relation holds 

CJ'C 1 = J, (Bll) 
therefore the correlation function can be reduced to 

- 2Tr (JGJG) . (B12) 

Introducing a vector potential which is coupled to the current vertex, we needs to evaluate 

- Tr (J K AG(J K A + Ci^AfC^G) , (B13) 

where K = C for the charge current vertex and K = S for the spin current vertex. In the 
case of charge (J K = J c ), with A = (A°s + v4 1 si)r 3 and in the case of spin (J K = J s ), 
with A = (A°so + A 1 s\)a Zl the following relation holds 

C(J X A)'C" = J K A. (B14) 

Applying to (|B13|) the second derivative with respect to A and to A 1 , in the following 
combination 



d 2 



dA°< 

we obtain in both the cases 



d 



2 



A=0 dA u 



(B15) 

A=0 



- (Tr(JG + JG + ) - Tr(JG-JG + )) oc a, (B16) 

which is the charge (if J = J c r 3 ) or the spin (if J = J s a z ) conductivity. Instead of using 
()B15|) . exploiting the d-wave symmetry, it is straightforward to show that Eq. (jB16|) comes 
easily from Eq. (jB13|) taking simply A = (s +i Si)r 3 for charge and A = (s +i S\)a z for spin. 

Let us consider in what follows only the spin case, with J K = J s , Eq. (jB9|) . since in 
the c?-wave superconductors only spin conductivity is conserved. If we now use the gauges 
A = (A°s + A 1 si)r 3 and A = A°s 3 a z + A 1 s 2 we obtain 

C(J RR/ A) t C t = -i(R - R')(t RR , - iA RR ,r 2Sl )A, (B17) 

therefore the correlation function (jB13|) becomes (5 — R— R) 

2Tr(5(t + iAT 2Sl )AG5tAG). (B18) 
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Taking advantage of c?-wave symmetry and using the relations 



TqS 2 GtqS 2 = T 3 SiGt 3 S!, (B19) 

t s 3 Gt s 3 = t 3 s Gt 3 s , (B20) 
which come from the following structure of the Green's function at fixed disorder 

G = [(t , t 3 ) <g> s + (r 2 , n) <g> si] <g> (71, 72) + «*S[(r 2 , n) <g> s 2 ® 73 + (t 3 , r ) <g> s 3 <g) 70)], 

(the terms T3S0, 7iSi, Tis 2 , T3S3 disappears if time reversal symmetry is preserved), we can 
find that Eq. (IBToT) on Eq. (IBT81) gives 

Tr (8tr 3 s G5tr 3 s G) - Tr (5 tr 3Sl G 8 tr 3Sl G) . (B21) 

Reminding that J c = —i8t and G~ = s±G + si, we can recognize Eq. (jB21j) as the Kubo 
formula for charge conductivity. 

We have shown that the gauges A = (A°s + A 1 si) r 3 and A = (A°s 3 cr z + A 1 s 2 ) r 
are equivalent to generate charge conductivity. The same calculation can be done to 
prove the equivalence of A = (A°so + A^si) r a z and A = (A°s 3 + A 1 s 2 a z ) T3 to produce 
spin conductivity. At the end, the possible choices for the vector potentials, taking into 
account also the replica space through a symmetric tensor, \ s = -^=(5 i5&2 + <5a2<5fci), or an 
antisymmetric one, A a = -W(<5 a ]A 2 — <5 a 2<5&i), are collected in the following table. 





Charge 


Spin 


A 


A s r 3 s cro 


A a r 3 s a z 


A s t o s 3 0- z 


A a To S3 CT 


A s t s a z 


AaToS (Jo 


AsT 3 S 3 (Jo 


A a t 3 s 3 a z 


A 1 


K t 3 Sl cr 


A a r 3 s l °z 


A s T" s 2 cr 


A a TO s 2 ° z 


A s tq Si a z 


A a T sio- 


A s t 3 s 2 a z 


A a T 3 S 2 CT 



APPENDIX C: THE ACTION WITH VECTOR POTENTIALS 
1. Charge vector potential 

— # 

Let us introduce a slow varying vector potential A oc r 3 , taking A and Ai from the table 
above. The fields change in this way 

CR _^ e if Jo* Ci? „ e i|AB ^ (C1) 

4 _^ J Re -iif R A R ,dR> „ c t ?e -iflfl > (C2 ) 
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since A is a slow varying function, nearly constant on the lattice length. In the Hamiltonian 
this transformation is equivalent to changing the hopping term in this way 

W — - t RR , (l - £f. A -^(5- i) 2 ) , (C3) 

where 5 = K — R 

The interaction term, before the parametrization by A, is unaffected by this gauge trans- 
formation and so 

A^— A^. (C4) 

The expressions in Eq. ()C3|) and Eq. ()C4|) represent the U(l) symmetry breaking. Now, 
returning to Eq. (JSTJl . we should consider in TrH rr ,T r i the transformed hopping term 



T R t RR , (i-itj.A-^s.Aryi, 



W (l + TrI VT+ + l -T R {5- V) 2 T^ - T R t RR , [i e -5- A+ ^(5- I) 2 ) T Rl = 

W + t RR 'S- T R VT R + \Y, 6 ^ T ^ T R ~ ^RtRRfS- ATI, ~ ^ R t RR! (8. A)%,. 

Besides the standard terms (the second and the third above) which bring to the non-linear 
cr-model as seen before, defining t = tiji + £272, A = A171 + A272 and G = g + i g2 ^ s2 a 3 
with g{k) = — g2 l S 2 [(^1 — ^A!r 2 Si)7i + (t 2 — iA 2 r2Si)7 2 ], we have the following additional 
terms 



1. i-Tr { Gt8 ■ A 



2. i-Tr {gTU ■ A5 ■ VT+ 

3. ^Tr (GTt(5-A) 2 T^, 

4. ^Tr (gTU- AT^GTtd ■ AT^ , 

5. i^Tr (gU- TVT^GTU- AT^. 
The first term is zero and using: 

i) Q = ftS, 3 T, ii) g = 1 -{G + + G-), iii) -^-^ = i(G* - G-), 
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vi) Eq. 



v) the c?-wave symmetry, implying that odd terms in A are zero under mo- 



mentum integration and finally vi) the relation (tiAi + A 2 t 2 ) 2 = (Af + A|)(if + t 2 ,), due to 
tiA 2 = t 2 Ai, we obtain, summing all terms and multiplying them for — ~ (the coefficient in 
front of Tr \n(e — H° + iQ)), the following additional contribution, depending on the vector 
potential, in the action 



S(A) 



TV 



32E 2 



:<J r Tr 



VQ + i-[Q, A}) ( VQ t - i-[A, Qt] ) _ (VQVg 



with 



a. 



a 



Y 2 

nrV ^ 



(V fe A A 



2\2 



(C5) 



where a is given by Eq. 



a. Bare charge conductivity 

Let us suppose to have A = A°sq + A 1 s\, to recover the Kubo formula we have to evaluate 

' d 2 lnZ d 2 lnZ s 



dA ' 



dA l2 



(C6) 



A=0 



with 



Z{A) = J DQe 



■So-S(A) 



Since S(A = 0) = 0, the generating function at zero vector potential is again Z(A = 0) = 
Zq = J DQe~ So and, introducing the notation (...)o to denote the quantum average with 
weight e~ So , we have 



d 2 lnZ 



dA a ' 



IdS{A) 



A=0 



\ dA a 



A=0l 



d 2 S{A) 



dA c 



A=0l 



f dS(A) 
V dA a 



A=0 



(C7) 



The first term is zero, the second is the average of the operator 



d 2 S(A) 



dA a ' 



A=0 



e 2 7r 
16c 2 E 2 



a c Tr ([Q(R),T 3 s a ][T 3 s a ,Q(Ry 



while the third is the average of the square value of the following operator 



dS(A) 



dA c 



i- (^Tr{VQ{R)Q{R) ] r 3 s L 



A=0 



At the saddle point, 



c V8S 2 



Q(R) = Q sp = Ss 3 



(C9) 



57 



Eq. (|C9|) is zero and the action is simply 



12- 



The bare conductivity is given by applying Eq. (jC6|) which yields 



d 2 S(A) 



d 2 S(A) 



A=0 



OA 12 



8vre z 



-Or 



(CIO) 



A=0 



where a c is given by Eq. (|C5|) . more explicitly, using Eq. (|8*3|). 



Or 



1 ui 



(£ 2 + £ 2 ) 2 J 7T 2 V 2 ' 



(Cll) 



At the Born level the c harg e conductivity that we found is in perfect agreement with the 
diagrammatic approach 36] . 



2. Spin vector potential 

Let us now introduce a vector potential of this kind A oc t§o z . The fields change as 
follows 

c m -^e^ M c Rh (C12) 
cju— e-*^ Cju , (C13) 

The parameters of the Hamiltonian becomes 

- in* (l - ^- 1- \{5 • I) 2 ) , 

A^A^i(l + e-^'-*>) * A^, (l - A- \(8. A) 2 ) . 

Although the parameters change differently, the two expansions to second order in A are 
equal and the additional term to the action due to the spin vector potential is 



VQ + 1[Q,A]J (W -|[AQ f ]) -(VQVQt) 



(C14) 



where now the bare spin conductivity is a s = a, exactly the stiffness of spin fluctuations, 
given by Eq. JEBD- 
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APPENDIX D: USEFUL FORMULA 



In the presence of chiral symmetry the gaussian propagator becomes 

{w; b Sinm {k)w; d Sirq {-k)) = i(i - (-) p A n A m ) 
:-) p (±)^ m (^)(c^(-) p [±]<e(-) i C-,^ 9 (-) p (-)i±]e 9 ^ 



ad cbc 
q^m—r 



+n nr (k)5 p3 5 S skoK b J C r d q \ , (Dl) 

where p = 0,3, depending on the 7-components, S = S,T for singlet or triplet component, 
(±) are related to real or imaginary matrix elements of Wo (the 70-component), listed 
in Section IXI Al [±] for symmetric or antisymmetric matrix elements (always referring to 
Wo, the sign (— ) p takes care of the sign differences, when they occur, between the two 
components in sublattice space, Wo and W3), (— ) % the sign that W acquires changing the 
signs of Matsubara frequencies (this occurs only for modes proportional to T\ and r 3 ), and 
finally 

D °^ k) = 4^ D k * + l-e m \ ' * h < D2 > 

Ilk 2 

U n , r (k) = D 3 nn (k) 2v{Dk2 - 2|e ry , in the replica limit. (D4) 
Using the charge conjugation of the unitary transformation U 

nayU'nay = U\ (D5) 

we can derive the following relation 

Tr (u^naU^r^) = V [±] Tr (u^naU^r^) , (D6) 

with 

— in p-h singlet channel, (I = 0, 3, o = <j ), 
(— ) l in p-h triplet channel, (I — 0, 3, a — a), 
— in p-p Cooper channel, (/ = 1, 2, o = a ), 

where [±] refers to the symmetric or antisymmetric Tj-Oj-component of W, with i = 0,1, 2, 3, 
j = 0,x,y,z (the same sign, [±], which appear in the corresponding gaussian propagator). 
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For the following derivative operator A = VUW, whose charge conjugation condition is 

T^OyA^xOy = -71A71, (D7) 

where 71 is the first Pauli matrix on sublattice space, we can derive the following property, 
under the trace, 

Tr (T Wj ,^ 7 X) = -(-n-) q [±U±h'Tr (A^a^a^) , (D8) 

useful to evaluate (S 2 Sl nt ) and I^S^S^), where 7 g and 7 P can be 70 and 73, matrices in 
sublattice space, while Uj and a,y are identities or Pauli matrices in spin space. 

Defining the quaternions fj = r , iTi, ir 2 , ir 3 and o j = cr ,i<j x ,i<jy,icr z , we have also this 
sum rule 

[±]ijTr (Mf t a 3 ) Tr (Nf t aA = —ATr(MN). (D9) 

where M and iV are generic 4x4 matrices. 

Applying the relations seen above, we can obtain, for instance, the generic expression for 
the mean value of S] nt in the presence of chiral symmetry 

(SL) = ~ E ^ qD tm 2 (l " H kl A mi A m2 ) 6( ni Tn 2 ±n 3 - n 4 ) 

where we have dropped, for simplicity, the momentum dependences and where Tr means 
trace over all degrees of freedom except in sublattice space over which we shall trace by tr. 
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